


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1991-09 


Performance enhancement of the NPS 
Transient Electromagnetic Scattering Laboratory 


Bresanl, Aldo E. 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/28445 
Copyright is reserved by the copyright owner 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 






































































































































































































































HGg@. 96 Ud. 0. CoM ae se 
Aw oh 9 O51 ST BAO a Gee, we. Scan mone cane fares bey 
Sayre Sey Pole aL ALS or, Lee RC end ateisha cea a! 
a oo abl Bei ata faty His SEQ ata Be ROr Dad «34 qtr egs ax Aghi nas. ce¥ brs soap east acme i 3 Ta Taenna 
, oe SS oN et 88 8 Re acy Sere et OEE dees ® ce rane ete maga es 2 Any eee ile 0A 9 ta aiviareet: ae oes. “Alesina eg pect bieataaaa: ioe 
bs tral bed 4 Rae “ 4 vel " a, ae qh aette 4 1 . ary midgigathovan css ; ina gaat 
i 2 ee a) ae sus Hy tet og 10 ere marie PY ee ee ee ttn 6 a eee 2 Vale be 21a Ke -ASA TG tay Picea a ppc pale A oho tiatetaen 
SO MC ge: 5 pager cope ate, oh Fo o\ay 0 thay 1 BS Wx tan settota anes Aaosiet kag liaties eh ‘es isis eae che yaaa a at sar leone 
; . Tans a es oa oie 0.3 TF bg dda Monit 4 O52 fe Pies Px pares RES aasives Viet G08, om PAsad st NE Koen efi nd seisnveaah oo Meike ant faee ere 
are ae oy Om a owes set 8 omy wae yp 7 » ee ae we! Peery | e'va ip BA ven nia a hive a ois fete na At. Poi d oy fey eienesea oe) : : a indie k Roa ote 
Hes ‘ene . 1 ‘ ‘ i "pe mu, Ay “ig é Pe va “en Hot on bdewea. ee ae 8 ty tds Ni es 1a'e wad pat un {01 0e Ste pelt Wa htoudor 4 ns sin nviniaionyt oe 
ime te ‘ parade Se LE YE " regtttea 4s oe ‘ese SEN bIA CASiIES Va cate tert ty wa fatwagiptay katy he pie. pains tde, Spite aod hight ree feria | (eae) oe papa Moots mala Seas 
‘ 1 . le ‘ : ns Van, 84d tee ] Vale thelry 4 2. 24640 O10 nd ! ty eens Te 
‘ 7 a open fb 8 epee 8 oe See a Ree gts te. tA fay *ytyhs sf a 1 Vass Sea Sinha cai a rylthe 0% ou 1g Far in aes ies giteatinspr ayaa 
. of ee is a a? *e,1 4 ““ E Boars at 4 re ‘y ies ante ie Pre | 4 Lice Oe ae cs plats 5 Nath bata pairings “ fis ee thee sta 4 eaten hg 
4 pe hoo Here 1 ae 1 1 . ' ry afe AA + ane ty tyee tat rela wt ay is a aa ane TWA os tit * rea hh rh Ph whee ti sera wih mi ae ats ah ey Ghires: 
° . a 4 ote? . ro bore Se ee oe ee ee | Hae Pes deo Phel aeper td He ii ano aha a 
4 . a a cy Fe ' ‘ S a he he ea e t 9 fe ? 4 the we wd it 45° Im, 9,80 ete gab eal tale erecta os 
t . . oe ee | : by 9.9, Rowe os me oma mS 6 AB Rt TPP onendatueaae Drag RB ens a 908 aM 
' ' ey tf on 8 ' wad vee 2 eT ee, se a” » it’ O54 249 hel eeu eae pasate Be 
; i : > eae ‘ab teas Ie N ia: Oh ee a Le See a {hd toes ee Oe eS ee Ota, U8. pea ALG, Ol8sBI wate! PLE CD neat (he hema: bas 1 
. . net 4 ar fe . 4 1 4e@ Coeds @ - Sigh aan orsHA 1 eo ee taeeiaale 4. Se Wiced saan GAs fae ».8 A ar Hein paanadongs Mecca teats. yamene fegacale Auer ve ahd aieey ment i 
’ ' ' b mer NW nL ela Gio n iN a OG 08 Ue Bart! Ae ee coe ae a ine 59 Saher yy Se qtsaahe aSsAs ‘ada Stats ye oats aa he orp! re depp s pAneasia a venciuaek aay wae rhecawa 
a ‘4 be anal Ce isl Oni aie . ‘ oe OF, ‘ ' ay ee | jet Agiat 9 = ale ty ., dike Hate. Mote Sedae »' it; 
i. 9 ee Gt o*, . err wi * A Dae ere Tate tae ae Ve/erasae it eae 98 faresdgace oe CS ake pratktnch-wsd weed a Staion 
' ors tb oyay se 4 peas teu ] tp wt fn ie vy ui . Pata? e WAP hg ee ae a = OBE te e' (Wes had he rir pres rm _ z waa aT ae 
; 7 x Se eetue aC Wd 5 Pn 800g My sae aus ft! we ' = bien doaies. Mise 5? sma acca & vee aan i trociger! ee ess 
1 o of *3 Mb Fi pore " - "oat, orp sey eV A 100 Gt ty? veda 98 sa aeaieee wiht ead ‘Sruadietehy BOs Lage oF SEAT ede edie Pie eon Tin Creme te wicy ty dk tere gt “. \ Say ite 
’ oS anes Jn UP ere ak iy eee ata ma 1 Ri Sas ta "Aiea Aas ON; ee Poe rery fa aL bs A Andsduiqar age aha Mh6ye, Ag dge Ged WiAshS Alas athe bak eee 
. . * ge meataitee . or. ve eat? a! te! oa i on) Tint @. a6 th ated Phe ers Stee tan a . pea bos stad AW ee brea oh Pe hades oo Recon “rie narefasteg 
eC et hee ee ee ee ed toa ewan ¢ Sava ie had, ot nt ata, nytt Pen tients re chien tad GPa ye pa casctse es ee ceeen a rate 
0 % ‘ oft bd ' + a4 ° weaned " a S'S rq ryt pees 0,9.) Ste, et 1 ater } Ab sm 8. es ale FA retell as UR se eeaidess cantatas % i cea sitar 
P age oh 8 ‘ ‘ se, 7 na f ih ° * pare e Rasta? oe : or ret rout ates Sy Adee ae: 5H é ' ay RR Saat ; pi a a 
: Pe os 310 i , peau : OEE the ‘ J % wird eo Nar atte 47 Al dnl 92190 Weariinn dg? J s a fe 
: 3 2 ie ne ee ee en - ee men x é a gaa “ara a : aie ie Sesteal beg Pe ee pes eyty itis ‘Nout 5S Iaracte zhane S94 eer +7) peeedant 
; m4 aie hay 4 Sm Wa) ease are eh isd! Shokeiew t Bite *honbee, ys re hls ob lb i tat hye A oh plan di hae Nort nee ott. s wet: A. armel 
r : * 1%. le . Ar p * 
eee eave wees we late Abeta, ‘Ped # Tt Hespg eaten A eee Rieti HEEL ER ae ah Tion aa erp ta iecaea 
ee c A peadaigas . : 2 i : denye 
oy ; nae Peeters ca ia oe te 4 re bret “Nas ai ved: ee Rt '$Aata Pre ip its hs Saw ateostenttien Das ecteike Wimiate aC hae ada alate ah a bo AC ne a hy, 
Ta See Yee a OPT) eae ate nee raged re t.,) “Nie Pensa ted begets hint as [Med eth 39 vet wae by pdPa! eatae 5 aerate z sme cate ae ois ecg Be ip nds ssi naeaen ans asad 
Mute eng Melaes Wy fae eh 2 Tegel a areata hort she f ela) Ta Bay ane eee Rate Add a "ha bp" 3 aoe ‘e, 18 ost vee mie ahs ee 
ne ir RT Leer Hee Nga ties ni Ud doh sane TRe crane tlinas isJate ls aiess Wagitwra ifranare® 4 gary ihe, ca Niedcaqatetn a8sF-A jaiccanblhaa: enatack 
' ' le : at 1 Sieecd 1 aatgh in {044 Liat y® 16 £45) . ae ' AR rate ee pe a a 
oh Pr, Se NO TUR or * one a st FY Rog ROL : as{Ale na cag, 01314: 4 rit wage ee {esbabs ts ea tad yes Ras coed sesh Suen a tnt vty nee forte 
F "s : a i OC en | the af aay mat 4 oe Pe eee ee - . ¢ “ ote 7 sea J a ‘ 
| ys ne 5 hi c : w "he é ind se bihatea ye us es ah aug ales pha be vas atic sotad weaca toshthed seas ‘. mae 24gey: Nicagin sta bait te Rerpaes yy iad a9 a0 
- fx ibe ere 40d 0 8 oda At he a rs osha py ie te we ytok ofs.a, sits Acai RCS ary A deainiise i ' rioting Pri crieoee : ce oaeabanaee ors 
migns eam ee Nd arts ‘4 ‘ 4 ves va. eray ary ate rise? Cie fea asgabivan 4 Shaltst ers ph ee a tet 
Of eet teh we ade ne 4 aon hy set tap. pints wah athens ean” (Make dry she: bar aires Ses “aes is 
. ; rr teage dyad as oe $ & dens Lo a! te + vqrcashzieelegs Jay Bee tard eum. ies Pe petenay ets pee aa ea sir salamat 
y i ; aL =e J pops 
bh ore q' 4 ie mee ee tg Sete 1] a, en tee: sm gisee nnd i mass wie Pepe rts elt esas ea ay : heh 5 
. erase} ra A ae tiAdetig & ue * Heh onda chien AWA Vent: cotete satis neh eisangachete aly Sagar 
BS Ye iwvaiccat a: aon ? 
‘ower ane en, A tut ae eat aiesigs oly Oh oe: 4idty 7 Bs ves es aten aya iaevacagmras ye pa! ReneS 
v4 cr. we Aigval un i A “re ite at on ah h n Pe eae areerips 
Ve et Gaal 4 ry 09% 5 < 4 oT navieg Enns 4e ee BERG owt, ie "aie fa Ci ano amt ae : i 
Ak Pubinars heared ved * wernt cg h her Binet wht» 8.) ae Mogg En ong hy ae res Bs MA Ys i aed co reale ye 1 GS Alin ma catdan hi asngn ss wel Seige sat 
Yeeaven piv is ‘ MP are yo fe as uf . ae an tal ” aU ced EF fi Syhiends neh By | al ae ie ale gage! siete si pee ete ee = 
. ' oat biv\eha ’ igs Ay } 
paren Taba Rat) 28S aah ¢ is eee ir arith An .80% age" Mea 5 epee hee! ! ar twats Oi cogeypeese : 2: re “dahl creer wee 
. 4 a wn Pia fame sone te Ai tat AS! a piles, crsal oc Aute 496 ¢¢ 2: tne tw pty ae PMePOeeT isan as . a - ee bees oe te ere, ae 6 
"ee Belo tet ne yn) ah een Y B.S rhe, 81 ce sar. wed WS i Serre atts Bree Hh ee ats bade aaah rs ae M fe fing pt 
y oa eo sett 6 +a '¢ y n9t "Sp ’ FAL ak on, tae 1 8 Wigs Sane Saat : ep tyra gy 
ae 2 iy ape ; t of + na ‘2 : Sto pas Stee wane Ht Hap arte id AER IL eoctg th & Sean 
‘ 











ee »:* o28,® 
eRe A wads att «ob 








, Dp mLaT ae ELA and By a dam ob Of 0 s/o AL Gehan 10%5 RAL IS Bots 
fs Meare tol cM lay ding jal yge dine “sate 4 ‘ Oe 9a td) aitke WA wt Boddy 
' ae 4% Pogo mi beth pee 19 esos bart Ere aap ary oro See ge me rote Pte slo ae asked as 

° *) fan ' be vine ° lv Oe Py Pe ; Seng bef iat kes 
oe gee ct; Maen depuis @ aceahedg *. Phi eat a vet ricer a hae swale ae stthi ghar ath’ Ds. atAt ehotais rg tow 
ad's oe Sang >’ RY aoa rh “ jas hc ah haa vyat E 23: . ri ' 

Veeafe "+ . ‘ 






















































































































































°% ae eR pinl be, Oh ge ah 
Sp 5 iget des mas ve rat} $1 nemy Sere thie. gl af ated ate & at 
t ‘ per ity * wall a asses sate’ #y* Ageia HY pease to i Tae zit iageite se Aphaee Soee ibsiat, bath ladbig to eaekiet saan ts as 
ine & ae te Ps Ce . oo ay ¢ 5.0 {4° fetes defo gcd a ecauuty pee it ai pat ees + Pot 
Cente te hat Ons Sot 4 abe Feat ae cabo v $:' td iy 13 onal t Paice pl a Tas. 3 vd a “° io <> Me Yrnatine ait bt ie ery Ace a arth iS gens MB 4 
de 3 ye » * as a) gly De pet ak ob 4 B puteon “gf, | rips ai fo Harr er i a X] ailys orien ms Wrest Ba a7, : mney ry) tes oe Rises met age 
‘ ? “Aes e gt i a tat 1 tee at Lbs wy opty’ sk tee Mah) cad it Rants *%. ih nt) aa? bate? sous Peel esty gee MAS, Rekenakere van mee 
i v"d erate . wail Vy vee 140 29" LS. Re Be Pebonet 1e¢5 FG carurs q ‘4 of Maye a ‘ é eee ws Wushe & 
ny oe v StS i. oh ' See by n ‘ ss a Aad ryt ata e.d telat vis > 4 city poet rh Stones i. SEE geal? a atce ek BARU ont : 
re ane ‘ ingheL aigd  A thacy eee ; ary te : Nasi Seempbed alec 
| aig nett ata lig apis hans etiaebe aga satan at i Pere ct ah ea RR a eee 
3 Ay i fee oe & a ae © StL soe mie : tet oy 4 ei 4 " ara Regen? Sie bed LE eR Ee nbigntardoniat 155 SNe geliets iit irae ; ee ms eal be emir 
é “2 1 08 © 0 Shyer i tte a ets ach aoe recy? tela, " aah, Ag tant ° C yee NG & a Leu ee sia . ehine vis aaa ‘bah penned aS a 
: Sa Nd . be BS petit a 7 Rely Sere en i 320% Layt' mele meing ApteA. og Se Saad y Reteeparaae iy 8 sae Neat at nis 6 Dan thy ying Sat 
eae : ve 1oohe Met HN muse ue at swore Pp ze Fok Auton rks ma eee gy , fant a 5 appellee oom Peat Hae ne Pesan bem, het na het Seer! 
wae L r ° ‘ ts t mm vue ‘ ' tity en te -' Call etd ; i as Ste HL i ts aay ETE Peinel's a hala Deets aria pe lagadtwacbeme 
4 \% 8 oe Ds ' 1 ob Ad trae vA. ree . a ba a at aur ae “¥ rh sitios” iannd> Row.’ i. Mage Pianctabbes mreu wide setagelotbontee fret ee e rgb eb to 
AS CO pele ro, mgatnede Be'e Uigi ta netaend ‘ yeh fackarnas lie 4 fog a ieF aay taretashes sagt Mice cutie Latte iia th a aha le 98 Ae 
. rab fee 9.2 8) faa nay A218? Oy Hse “6 HES ab AES oR badget” fae Wp Bt pd Pel rat PrP Tree Aas Perk re: yy ahem, 
. ce dene re Slag eo pA Des ae tds Cebyl FAD ac 8 FCA, A D a THe erie Tey te ety 5 
: ou: : papi: <5 P Se qed Beef, ne ites pie % as ie *F ae eine Meni Rela Oia isle ean uate xe PLCs mya teneee 
Pesta MANNE O ME oi eis nicele tine. f5 8 3 ie are ren erat? Phe: a adeocct ta 7.5 rere 4 eas Ya Seee = tet ath by 
: B ere v ors reer 6 rie arse Ts oa rie eee Rica Hi dit A ME iy Hor y oPenta seks fy 3 Rie eae Pee De etek org bea taeaN cake 
, Per ry a o 4 af-ea al b ® 7. he - 1a} et : 
‘ BM ner) ees Be tem, ie 8h! he, Bn vistas sr ot hae ester re ses ava ed 952K crepe on REET LOO 
. oe ; Pe abe boyy EH 2 8° Ve o "i af - r a 4 
‘ Cero aet a gate Cathie vpet ow v4 my ” atts hate tata ears ath “hi v8, HH sate nagte Vege’, MERIS 0 jy wo Ta ak MG Ooty Bete 
See aM og Morelia es SE oe at Ait 82 athe ao CoS? veka OLE 4 Bae YP ere we) wena hid ieee ae <r 
sae : ANKE Cys Asta komt J dyge eaber YOM rag y ately a Agl thine a's ai al baa “ a y eatin nwioe wid af cee ot 
. yy a4 Me g® os eae wot ug 78 re ag gt | “py at’ as, Fipaidtedt y Rat sida J pirat vy 
Bre stot Tistasgh, H le * iy adhe ‘ie ata: wanes tt att asi ming i « Chetan 
. ’ . %. , s° a. pater sts é "5 ha, at q: 
: r dnd ai he ss hs Pha = 5 H me as 
1a? a oh . hs Bas V seu KM ty 
v bee, le He . 7,4 ‘ 
a Sreeemch es Lad tid Seba is Wel © ae wants cing wei re heen 
re 
‘ Be . 


ius hgh Rdg Male PS Rofo y Maly We 4 Me iia be, 
1B OE vobehg hte Sab teagetye mw bts fotebotas 
gages Gighonasken ‘ 






Sida, Aetings 
ta genee ie SERRATE 
tae lg Wha kota? 
ne Seeee nt oe fea ehnat tte 
Tn oh wth Seo Badal! Lt 
Eine Ra Sat "ony “% © we Aas afte 






















* 
°a A he (Sas be 


































































































































































































































































































































































































































































































































































































































































































tet sistiie eet tens ints 
i ies ris BT Pp rotor 
Sh, ' beak: Ae ir 
‘ o Fee Nae fhe EER wd Ske |, 
‘ i 16 ; yo Y 5 A Mien 5, Miep Lae aA pe 
De eis iy ce at Pare Be Foe atile ree ut aides rl Ota eer bes Ete} 
Re ia Lt ere) ae Sens nite | ray 1% oF atin, pes aL, rest Bees ce nat ine 
Pe ek ep eimenad yee te a Oey oe hee: “hes . ee Pere: entree b Ale ae saat 
Bie alee ty" ape : ‘wt ; ates Bi oniies wide bars te theta 3 sae 
a) tke ae . Ak oeracan ' is . : Sdaee VGGRS Ea). Mose : £5926 Oo Pet Oe is as Hons 
‘ , cies mee fit ota dow ee Beane ca$ yas tepeetn da ALE ge? Nee finbe Be oe pet teyage: hn ae 
Ue Da ces ee ' ratetod tot Ome Bo surn wae Urs ; a spigaer 4h H ee 8 es phd eee Sea it eee . 
Bie sites pe fe nee cree Set recta oft! eft ve . t soefte is fed 
eat Mt uote Hho af 4 ; etsy SiN G58 gal Tea SPAN GL Su, eas i 
rie : tilea hu Watney alas, sakfeatt duke it acon beh en ae 
tee ? ATH Tory" EE gs os sg 10 oe Mae fad NSP RAG Ve 
‘y et 4 fw mii 4 pie Squaw War Sik Bioy” Pe acter ee 
Le . : _ Peer ate con eet 
a ied 
oe, . S Heng fee 
uer mare te oat sete AS ike ditt 
4" ihe . G4 asa sah ect 
sears eae aH nS eeeh yl: BENG Siatice si * woh BF 08 “al He ABTS Abn RS : 
‘3 ( es eae tides Mi wet yeas vise EATS s hg hs) earn ize oe nsbatee Sag? FE 
' events QeleSiE + vei braids whe He nares ge 407 fre eh =" bre a es 
Coe ea as @3e : fe ‘é ae 4.8 eee ee 
5 : : ike é “i peewee he = Oe EB bt 
a ere r ty. $3405 oa. RBM m hh 
ie Rees fot 3 4 3 
e 4 
+ ae <f 
bie ial" 3 3335 i 
. s wh 
e « Us s é 
ae wo rhe f) trail 
ee ee ‘ ae 
Vodtytet ya ae 4 % ° 8, 
‘ a. : gt “ 
iJaeole 32 é sie “s 
ary af eft rd ite 
sie a es ' 
, tatoo Sofas 23 
oy 18 ey “$ 7? ; fas 
as ie afi ie 
7 * "e« ow Ne 
ar; pace gaé a $28 Oe e 
ya teat 
gule gk A FH, Pass? 
ese —s oF Iyatf a 
wt.t sahas eat : hs 
e a2 23) 
‘ae 7 Paes i ; 
sae ep Sled is Mae oe 
~ ro . 
PCAC) E is $s ee * 
. 2) aN, 
ee ot &: wits "ou oe errthh 
atti tee Ay Py qt: delta soot as HAN 
es : ie eed tf Pry Ty ion, 
eee. ete pt al ts t Py ats 4 
. 4 oe. M . vq wena a 
reg ot be ; oF 00" te ry eet 
ie ; : 1a in ek hae Ne gu ae 
A : ig ‘i ‘ : sta nokia ie eee das f sate 8 page aN a3 xe ‘ oa rate) ry 
il Leet i of et bt oa + oe rah 4 oe 
*.% ot ot e : ¢ 4s in 0 e0,-geum x ‘3 Ly i Ah ic ta ses . 
Teer, aie. AL “fils nh 4 ene ie, Tete e 
pote 2 4) Fret 1 Ae Pa regteee Pay “3 pionee see oe Sitges TY REA eee 
te, fe te tee See 1, t F ars . "rAtgaw. et whee gianna ce) 
ys . owe 24 sete Mt i teprka tao sai hase 
a. tae “ 3 vr 5 4g 
4 a ’ , pie fey tae ' ° F st Rie r% 5 a7 3h AveNg erie 4 
Cop ee. re era > on ta +. yal ereety ade Os ee ae £§,, USeag eR? att fiyek tae pares 
© ‘hee bd eS ae are ee at ts ae ’ ‘ E | &3°, fe over ares Pama 
ores me ar | Vig Pat: 4, : 1 0°48 gg FF a re) nv) rei ic aie ae Po Peteh HK ates eres 
. e 4a of  w 7 ot? 4 < pater d o oF ASA in HRI tales 
. . ' Andee # ete esdrogrey, toe shoes ’ 5 i aes eh 
yA) Ser eee, yo - Me on on oh Pt Ss | ae i % ; ee * : aetet ‘ Farete, Ae r stances , Ee 
ees. eee ene Uamererica ya & ane mg gees? <b Ss ah a > oe ots es eye tea (he vat 
oot %, P ove eg j aE eed da lisge ct ys cae ayn meee by or ° rea a Lele a)y vee ri he best oe ft tipive ote ° ne (pee aa ove 
fa ' 4 “rat t eae os. on Vy a a § at * yey 4% fos . { He Leg WA Aine gel ge He at ¥ 2 
ate Se Se For igs bike too r) . ’ Da we air Fd ay! er ert aera vege take ae 4 thsireae iy ‘_ 
° t « . . = pte A ar ct oaepe as 
Cte I ae aeigeientii so eiie ie . ke PS hath te Ebest - tee see ugiverics meters’ gern 
or s es ; i 33 “ ‘6 . oF Ee Pee tee OC pi br tye uae fied wert 
5 & wo) 8) tee yee . WePeke re: oat © 
' ‘ta _ ; nt rely es an: pga ae of ore ig y 2% H ere hana Sena 
a $ : a) eigen eres oc rTny, ( Pre Bet heserespey bait beg iter e: 
' 8, * ee 2 tae beet : Les 4 ; aa . “ o teerls ra ore dt iy p bagi oany 
® rd Tt : ‘ A ok : ! mr. BFF 2198 poy 
Site Meares Ree CR Oa Peas ge ieee en ate. cee rie? * we an ta i orangt Sond 
OS ‘ ve pie ee a ree bee So tapas Meee Rae: pallens pe e413 a2 
fora ate V8 2 ey ' 1 fe Se ¢ itrgyr oa coaeee 
petit pe . om, te, ' ath ha tae Ase 
Boa Py ;" ‘ PTT eC ted he bd 
ems ee es we iri [sot onvhgen! otic ota 
; P daa : a <C 2 te, vac th 7 + aby wer Ay we te Ah oe ay 
ae os es , Werk. a Ae ints $PhE nehad: A ei hed bid ba rr) 
? , 1 
a te 4 a ve ’ ea BLN ae OH aty mete Ae sid ite & fore aeed. sheep gach 
ter ae ter uth Harte at vite peepee) Bs Han ris canue ag Pete 
ote 4 J . © 1 sf 
i aT . } ho ed ye ecek gre? {ize . ce F] tide 9 at eee Bogs 
& > ‘. 1 $4 ' Q¥. ae } rie sev nm hye Peed | 7 eta 
Ov uae H coe eas et wes : roar a oye ye NWP, bash aos AG or stb yo ie ashy mass ifey ie 
ete es it) a << we Ai " iM AyD n genes r Coed be oh ty & oe 14 t ey athe so jopee satay 
De Mae . e ae aT u idee we atts vie wee eae tae tered eit 
a : ; C Cit a ot faa! F " a5 Hy mts oe Gah kane teh re Litt soa wide PS cr iY *% op ay fy. isae 
one 7 i Es 4 Cay vate ne " ne far a ar “ tt eit “Sint An ae cial ers? i small! ad one nat castes meres df 
no eats ot dD aa ow - beriethy HOU Slt alot Tae Pantene kuti ert Med igen 
‘i . ae ‘ae press ess rete out le Me : p +e bs iA. es BYeAP reste! races cn pire Pei pret es Uhh eon Thy ae att A i ret ger, 
' ' : ORC a cn . “ 45 ms sv Hae fe 4 3 « SO Beigt. Yo 1 Venah. 6 eet 8) oe adhe we x at pimioet’ 
ae Ae: A ee parc a) : : ree: et Hy eh ie Ee Btn . Poa R oer *s at % Pata His Hart: vas ANA Rash oy mee Kon ap) See edt or a 2h 
Aare 1 bee i fo 7 Pe Bde eb igloo ote bolas te, vr ben bghe ! poate? ah * is om Heyeny m5! seveed ite ura hy dase Hin Jogi =peee Sian Gite ee or ra Tezeee ead 
: Pye oe ‘ha 4 Ge . sd ral doh . Mae ate ge Hab Sac fry a ‘dapeeseg a * ee vei avec als ved in sone natin gene ctuy Aer eek. case anh eye 
ae ee BRP Sas DTC rR on MRA os i eae dgeele ce gi ee ate a ee teeta Be Boe 
oiistittievers oy ast) weer ogee ped ee ee apie : HAD oP eine BEES aan Chen 
Ou Gee Oar ian |. NP I. SR ee aa Rae Ages! Re PG iis MALE ee Els tat ai 
‘ . . 7 ners, i . we seek { . *% \ 
eek . gh ore Pow 40 er! 4 of : { is Le err rae fh, iT. A% fener? it 
» 4 . + f 0” ne os. Witistes 
i oeowias ri ML pa AE eG “agtee Per) wee vt edits hts Uris pea ese BA v4 
oa stl oe Crs § Hay & N (Aa ssetha 
aie 7 6 er 1 ' Stare 5 9) et ad abe, te tester, t6 
> oe A aaa ee ot te pas! ‘on agi i lin dee ik ae Marais a ies . Reais eile 
; Pa ie re H . : r; UC ae OF: 2 vide wetea 1. i re Phd ala AS LEA pe: ie ont S 
rare eee ot. a Segre Me ATi "i a, irej'pe $48 Hee 484 os ni art 
Oe ie -. Hye ae j a tig i oy vie as navi gt 2% ied serine ae tf mI Pathe 
se *" a ' oor rt be ; if ROK’ do Saat ane ny bo fe ade ey Tate pais 3S Pants treetarabagatly © Ohta ie we as a ae wise 
. ; "1, ° or ‘ - 4% Sb pet ng the Py ee a ha oil tage i 7, ened f yi Weer ashes bok 
oh it te voharl, oie Se ON Ae eae ce A Taya ite 4 rei Aa, hes Aathagen s sieges Uae die hte Stgnes BSc tos ae futay we vies f 
- aes , : me con aa! i, states Ait wae ts 2 0 arte ite = P(e e vugey. cae eat JH ” sijierare me 4 i 
¢ peers Sie PRO 9 ether acthria oe oe “i + Mey S54. sar.ast 7 nl ee 
i <r sp tge o's Jan £ wt 
xe ACR ey See sh i OF PS LS 2% e Mets t A Re "yl mire js" 
ae A. eo pe : > topler ahh? Sh a we A at Hee Ary telaby’, At ry spe see 
A er Pees. See re Natit ek “tee Ce yo ey toe Gucan Oe 
Ms ¢ ue oi * r a | *, ‘ . Pises o 2, 4 Ant i Pete ae 
pat sie 8 : . é A some gt ail ee te pUrtulelarstost % + 1a is sh ytel rn cane 
nen Per ees ey ae Ne iS AY iS rooate 2 byasaee slut gas oes ie eat pce 
=m ce 8 - te 2 % « H 4 LW £ 1d 
a ER EC Restate aiporn in ENC aE URAL euitidatatnaet 
. ae ane ab aeolian’ rae rece ea “ setts Poaatinan I fia Sk Sinstesa: gee CPF Uh ste 4 Lae hats a ie : ve way" eters esey 
' ° ‘ tp 09 . bt ‘abetect ; + > TEST LPL PR ea oy ¥ ee bt ep ei -giergs es iy 
Cuore °° a *y'  aihihy a cst das “shale ol a HG uty fetes Menge ke awl des Be Eis us us teed fin ethyeai hg Mets Py re gael ne fatten! rit oeoues erg 
. A 1 t rT F a ’ . rp ‘ee On + o. oe raZ bas % AD jen Pv 4 yautee ys ae MaMa Pad he of i a ge “4 uy ve orien s, a 
a ee 7 a ; . ; . ‘ 5 es a "4 8 (a oe 3 shor bi J j bo : ; vans rans ’ ore eri? ER hs HN aN ae sneer Ut Fs aa RAN ioe 343 “7539 DR cit tte 
Cal H AD 5 ' Ch % F ~ 4 aie GV bape F On Be in 3 ody agst oe ay 
' Popat lie Jee ‘ 5 . ’ of * "a es Lae He rear Ge) Perna \> atte PIO Wey; iw Aeon R ESAS L seep 
oad ne : ape es ve : LW UR LET LU as suttbog SK the sleys PAE "4 any ING: Wie ate Be gt Sarasa. ob ge wLy 
a ade i Be ete Gea Te a eete cd oe tates uit Cares eth ot enaes BS ey 4 ok PA et te lieren tit May seats Y 4s) Sapte, eeues eh vals 
bees en 7 hae tte Aone yy ante nace ameraiem Be she Bosh sa Stay ‘ te aM. ges barton age Heats a Hea free avin ceteris pay Mtge RAST 
oe ° . oi i i e . ae o * 
: : ay aoe eres”. ( RCL Ht YS naka bos eon tthe md Gee ng ies nse cae ees a 
. Oe ce oe eae eee eke : ‘ Pr o ohte 5 a jhe any hg NM he ta Sines ae Pe sae fy nbs: besa 
: aise h 4e Rartceare p aetae ri ats 
; 2 "8 od . ey iP 8 we 3 ha.) VY ete 1h) Be bt eatys, ; srs ab wy ey at An rove Pld Sierovere ; Jet enw 
@ ‘ 7 Re \ fe ; ckbten ee % Linea ° be + Ra, apa eet cvil ey) Up et foe iis Tees, Vader Py ke 
% 7 5 G te ' Fee i = 4} . s re htla Scsdy & 10k Leal 
: z) Leer tery 2 O85 ' ‘ 
. * * bd - . a 7 8s 
ree ° . 2 











NAVAL POSTGRADUATE SCHOOL 


Monterey, California 


PERFORMANCE ENHANCEMENT OF THE NPS TRANSIENT 
ELECTROMAGNETIC SCATTERING LABORATORY 


by 


Aldo E. Bresani 


September 1991 


Thesis Advisor: Michael A. Morgan 





Approved for public release; distribution is unlimited 


1253944 
1253945 








- -“©¢ 


UNCLADSOIFIED 


gp RR 


SECURITY CLASSIFICATION OF THIS PAGE 


REPORT DOCUMENTATION PAGE ae 


1a. REPORT SECURITY CLASSIFICATION 1b RESTRICTIVE MARKINGS 
UNCLASSIFIED 


2a SECURITY CLASSIFICATION AUTHORITY 


2b. DECLASSIFICATION/DOWNGRADING SCHEDULE 


4 PERFORMING ORGANIZATION REPORT NUMBER(S) 












3 DISTRIBUTION: AVAILABILITY OF REPORT 


Approved fur public release; distribution is unlimited. 









5 MONITORING ORGANIZATION REPORT NUMBER(S) 










6b OFFICE SYMBOL 
(if applicable) 
EC 


7a NAME OF MONITORING ORGANIZATION 


Naval Postgraduate School 








6a. NAME OF PERFORMING ORGANIZATION 
Naval Postgraduate School 










6c. ADDRESS (City, State, and ZIP Code) 
Monterey, CA 93943-5000 


7b ADDRESS (City, State, and ZIP Code) 
Monterey, CA 93943-5000 


8a. NAME OF FUNDING/SPONSORING 8b. OFFICE SYMBOL 9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 
ORGANIZATION (if applicable) 





| 8c ADDRESS (City, State, and ZIP Code) 10 SOURCE OF FUNDING NUMBERS | 


Project NO Task NO Work Unit Accession 
Number 
11 TITLE (include Security Classification) 


PERFORMANCE ENHANCEMENT OF THE NPS TRANSIENT ELECTROMAGNETIC SCATTERING LABORATORY 


Program Element NO 





12. PERSONAL AUTHOR(S) Bresani, Aldo E. 





13a TYPE OF REPORT 13b TIME COVERED 14 DATE OF REPORT (year, month, day) 15. PAGE COUNT 
Master’s Thesis From To 1991 September 104 


16. SUPPLEMENTARY NOTATION 


The views expressed in this thesis are those of the author and do not reflect the offictal policy or position of the Department of Defense or the U.S. 
Government. 


17. COSATI CODES 18. SUBJECT TERMS (continue on reverse if necessary and identity by block number) 


| FIELD GROUP SUBGROUP Transient Scattering, Signature 














19. ABSTRACT (continue on reverse if necessary and identify by block number) 


This thesis describes the performance enhancement of the NPS Transient Electromagnetic Scattering Laboratory (TESL) accomplished by 
replacing the old H1P 8349A microwave preamplifier of the dual amplifier configuration with a new Avantek 13533 5-13 GHz amplifier and 
optimizing the delay line length for the 1-6 GHz amplifier. New Matlab software was developed to process the signals scattered from canonical 
and complex targets. This software includes a program to calculate the correct delay line length for either amplifier in future modifications. The 
updated TESL is shown to provide measurements yielding excellent agreement with theoretically predicted responses of canonical targets 
demonstrating a significant improvement of the signal to noise ratio as compared with the previous configuration. A target library was created to 
support research in radar target identification based on natural resonances. 


20. DISTRIBUTION/AVAILABILITY OF ABSTRACT 21. ABSTRACT SECURITY CLASSIFICATION 
UNCLASSIFIED/UNLIMITED 0 SAME AS REPORT 0 DIIC USERS UNCLASSIFIED 
22a. NAME OF RESPONSIBLE INDIVIDUAL 22b. TELEPHONE (Include Area code) 
Morgan, Michael A. 408-646-2081 EC/Mw 
DD FORM 1473, 84 MAR 83 APR edition may be used until exhausted SECURITY CLASSIFICATION OF THIS PAGE. 
All other editions are obsolete UNCLASSIFIED 





Approved for public release; distribution is unlimited. 


PERFORMANCE ENHANCEMENT OF THE NPS TRANSIENT ELECTROMAGNETIC 
SCATTERING LABORATORY 


by 


Aldo E. Bresani 
f 
Lieutenant, Peruvian Navy 
Peruvian Naval Academy, 1983 


Submitted in partial fulfillment 
of the requirements for the degrees of 


MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 
MASTER OF SCIENCE IN SYSTEMS ENGINEERING 
(ELECTRONIC WARFARE) 
from 
NAVAL POSTGRADUATE SCHOOL 

September a 


ABSTRACT 


This thesis describes the performance enhancement of the 
NPS Transient Electromagnetic Scattering Laboratory (TESL) 
accomplished by replacing the old HP 8349A microwave preampli- 
fier of the dual amplifier configuration with a new Avantek 
13533 5-13 GHz amplifier and optimizing the delay line length 
for the 1-6 GHz amplifier. New Matlab software was developed 
to process the signals scattered from canonical and complex 
targets. This software includes a program to calculate the 
correct delay line length for either amplifier in future 
modifications. The updated TESL is shown to provide measure- 
ments yielding excellent agreement with theoretically predict- 
ed responses of canonical targets demonstrating a significant 
improvement of the signal to noise ratio as compared with the 
previous configuration. A target library was created to 
support research in radar target identification based on 


natural resonances. 
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I. INTRODUCTION 


A. HISTORY OF THE TESL 

Frequency domain techniques have historically been most 
often employed for electromagnetic analysis. The band-limited 
impulse response of measured targets was usually obtained by 
sweeping frequencies over the entire band using a continuous 
wave source; time domain techniques were seldom used due to 
hardware limitations. When fast pulse sources and sampling 
oscilloscopes became available, time domain techniques gave Sn 
alternative way to obtain the transient scattering response of 
measured targets ina shorter time. 

The original Transient Electromagnetic Scattering 
Laboratory (TESL) was initiated in 1980 by Professor M.A. 
Morgan who has developed it through the years with his thesis 
students. Hammond (Ref. 1] worked with this 
laboratory, researching inverse scattering based upon 
synthesized ramp responses. The original TESL incorporated an 
outside ground-plane range and was limited to symmetric 
targets which could be mirror imaged. In 1983 the TESL became 
operational as an indoor facility located in Spanagel 535. 
Mariategui (Ref. 2} measured scattering by suspended targets 
in a shielded anechoic chamber, which eliminated the 


requirement of target mirror symmetry. This initial free-field 


range provided good transient scattering measurements. In 1985 
McDaniel [Ref. 3] showed good agreement between 
theoretical computations and measurements of the impulse 
responses for a sphere and a thin wire. 

In 1988 Sompaee [Ref. 4] incorporated new 
measurement and computational hardware and designed new 
software to automate high quality transient scattering 
measurements. The major hardware included a microcomputer, a 
new HP 54120T digital programmable oscilloscope (DPO), and a 
shielded anechoic chamber. Automatic transient scattering 
measurements were performed using an IEEE bus controller. He 
showed very good results comparing experimental measurements 
and computations for a sphere and a thin wire. An improvement 
of 10 dB over the Tektronix hardware previously used was 
achieved. 

Prior tO the current effort Walsh (Ref. 5] 
incorporated two parallel GaAs FET broadband amplifiers to 
extend the frequency range and improve the signal to noise 
ratio (SNR) of the TESL. An improvement of up to 8 dB in SNR 
was observed when compared with the single amplifier 


configuration. 


B. OVERVIEW 
The emphasis of this research was to investigate several 
improvements in the implementation of the dual amplifier 


configuration for the TESL. Hardware improvements involved 


replacing the old HP 8349A 2-20 GHz microwave preamplifier 
with a new Avantek AWT 13533 5-13 GHz amplifier and optimizing 
the delay line for the 1-6 GHz amplifier. New Matlab software 
was developed to make the process more user friendly and more 
efficient, allowing an in-deepth analysis of the measured and 
processed waveforms. The new software includes a program 
employing two different approaches: time domain and frequency 
domain deconvolution. In addition programs to calculate the 
SNR, noise performance and noise spectrum of the TESL have 
been developed and tested. Another Matlab program was created 
for future changes in the TESL in order to calculate the 
correct delay line length for each of the amplifiers. 

Validation of the new configuration was done by comparing 
experimental results to that theoretically predicted for 
Simple canonical targets. A library of high quality 
measurements for simple canonical and scale model military 
targets was created. This library will be used for testing new 
natural resonance extraction algorithms and resonance 
annihilation filters as applied to aspect invariant radar 
target identification. 

This thesis is broken into six Chapters. Chapter II 
describes the physical laboratory facilities and the changes 
incorporated through this work. 

The theory of transient scattering measurements, the 
mathematical model, and the transient response solution are 


explained in Chapter III. 


The description of the acquisition algorithm and 
measurement process with examples can be found in Chapter IV. 
Also in this chapter, one can find a description of Matlab 
software developed in this thesis. Experimental results 
compared with theoretical predictions are also displayed, as 
well as the SNR, noise performance and spectra. Moreover, the 
target library created is listed in detail in Chapter IV. 

Chapter V presents an overall review of the work and some 


suggested improvements for the TESL. 


II. LABORATORY DESCRIPTION AND MODIFICATIONS 


A. LABORATORY DESCRIPTION 

The anechoic chamber of the Naval Postgraduate School 
Transient Electromagnetic Scattering Laboratory is located in 
Spanagel Hall, room 535. It incorporates a metallically 
shielded anechoic chamber to avoid outside electromagnetic 
interference and to provide minimum reflections in order to 
Simulate free space. The dimensions are 6.10 meters long, 3.05 
meters height and 3.05 meters width. Targets are measured ae 
a low density styrofoam pedestal, located 2.18 meters from the 
antenna array panel. The interior (floor, side walls, and 
ceiling) is covered with longitudinal wedges of a special 
carbon impregnated foam absorbing material used to direct the 
energy towards the absorbing back wall. The interior of this 
chamber is illustrated in Figure 1. Further information 
describing the chamber can be found in References 2 and 3. 


The hardware components of the TESL are the following: 


HP54120T Digital Programmable Oscilloscope (DPO). 
- HP54120A Digitizing Oscilloscope Mainframe. 


- HP54121A Four Channel Test Set. 


HP-IB Controller Interface. 


AVANTEK SA83-2954 2-6 GHz amplifier. 


AVANTEK APT-12066 6-12 GHz amplifier. 
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Figure l. Interior of Anechoic Chamber 


- AVANTEK AWT-13533 5-13 GHz preamplifier. 

- Midisco inductive power splitter. 

= 141 mil Semirigid coax cabling. 

- Two cooling fans. 

- Two transmitting horn antennas. 

- One receiving horn antenna. 

- PC 80386 25 MHz with 80387 math coprocessor. 

Further details about target acquisition software and 
hardware can be found in References 4 and 5. The exterior 
hardwareseompenents are Shown in Figures 3 and 4. The NPS TESL 


general layout is shown in Figure 2. 
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Figure 2. (From Ref. 5) NPS TESL General Layout 
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Figure 3. Exterior Hardware Components 





Figure 4. New Amplifier Configuration 


B. MODIFICATIONS 

The most significant modifications done during this thesis 
to the TESL were the replacement of the HP 8349A microwave 
amplifier by an AVANTEK AWT-13533 amplifier, (which is a 
higher gain lower noise device than the Hewlett Packard), and 
the delay line for the 2-6 GHz amplifier. The amplifiers and 
antenna configuration are illustrated in Figure 5. 

Walsh [Ref. 5] reported deep nulls in the spectrum due 
to destructive interference between the waveforms in the area 
of spectral overlap. For this reason he empirically deduced 
the correct delay line for the low frequency amplifier Reins 
a tedious cut-and-try procedure. 

In this new effort initially only the upper band amplifier 
was changed and the minimum length coaxial line was set up for 
each of the amplifiers. Once the first measurements were taken 
and after signal processing, the spectrum did not show any 
nulls but the time waveform had an oscillation in front of the 
response. The analysis of this oscillation showed that there 
was not any problem related to destructive interference. It 
was thought that the problem was in the signal processing. 
Several trials, that will be explained in Chapter IV, were 
done with existing and new software in order to eliminate the 
undesired oscillation without any success. It was determined 
that there was not a software problem. A Matlab program was 
created in order to simulate the effect of the delay line in 


each one of the amplifiers. Two sets of measurements were 


10 


taken, each with just one amplifier in place and the other 
amplifier replaced by a dummy load. When the program was run, 
it was determined that the delay line had a large effect in 
the time domain deconvolution waveform, although the spectrum 
did not show any deep voids. 

Each measured waveform has 1024 points corresponding to a 
20 Nsec window, so it was determined that the best results 
were obtained when the measured waveform from the 2-6 GHz 
amplifier was delayed 30 time points. This corresponds to 0.56 
Nsec. The dielectric constant of the 141 mil semirigid coaxial 
cable is specified as 2.0006, [Ref 5]. 


The speed of light inside the coaxial cable is: 


8 
pee 3*10"m/s = 2.12 +108 m/s 
where c = speed of light in free space 
€. = dielectric constant of coaxial cable 


The extra delay line has length: 


e = vxAt = 12.0lcm 


The final delay line for the low frequency amplifier was 
25.5 cm, which yielded the best results. In Chapter IV typical 
scattered and measured deconvolution waveforms of different 
canonical and complex targets are shown, obtained using the 


minimum coaxial line and with the final extra delay. 
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IIIT. THEORY OF TRANSIENT SCATTERING MEASUREMENTS 


A. TESL SYSTEM REPRESENTATION 
The Transient Electromagnetic Scattering Laboratory 
system representation transfer function is illustrated in 


Figure 6. 
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Figure 6. (From Ref. 5) TESL System Representation 
The following table specifies the transfer functions used 


in the system representation and the mathematical model is 


explained in the next section. 
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TABLE 1. SYSTEM REPRESENTATION TRANSFER FUNCTIONS 


H, (£) 


H, (£) 
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[rules generator driving the Tx antennas 
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' Transmitting antennas response 


Receiving antenna response 


Direct couple of receiving and 







transmitting antennas 


Target response H,(f) 
Chamber response H. (f£) 
Interaction between target and chamber H., @f) 


Target measurement Hy, (£) 






Multiple scattering between target and 5 pes 3a) 


absorber for target measurement 


Calibration measurement Hye) 


Multiple scattering between calibration 


sphere and absorber for calibration 





measurement 


Total system noise N(f) 


Final signal which is sampled vee) 
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In the TESL, the pulse source is the step generator in the 
Time Domain Refectometry (TDR) system. This output is 
available from channel 1 of the HP54121A shown in Figure 3. 
This sub-assembly is also the coherent trigger for the 
sampling circuits of the DPO receiver and thus provides 
exceptional time coherence of the system. The risetime of the 
step generator is approximately 25 psec and it was operates at 
a 500 KHz repetition frequency. 

After the pulse is amplified, it drives the transmitting 
antennas which radiate into the anechoic chamber. The radiated 
field couples directly into the receiving antenna, en 
interacts with the scatterer and the chamber clutter. Multiple 
scattering between the target and the chamber clutter is 
indicated by the two-way arrow in Figure 6 on page 13. The 
total noise of the system is due to the thermal emission 
inside the chamber and the receiver noise in the sampling 


front end of the DPO. 


B. MATHEMATICAL MODEL 

In order to obtain the transfer function of the target, 
three independent measurements are required. Normally the 
first is done with a calibration target, which is usually a 15 
cm diameter sphere. This is called the calibration 
measurement. The second is with no target present and is 
called the background measurement. The third is with the 


target of interest present and is called the target 
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measurement. These three measurements can be represented by 
the transfer functions specified in Table 1. 


A calibration measurement is represented as: 


Y,(£)=X(£)H,(£) H, (£) (H,(£) +H, (£) +Hy,(£) +Hyge(£) } + N,(£) (3.1) 


A background measurement is represented as: 


Y,(£) =X (£)H, (£)H, (£) (H, (£) +H, (£)} + N,(f) (3.2) 


A target measurement is represented as: 


Y,(£)=X(£)H,(£)H, (£) (H,(£) +H, (£) +H, (£) +Hyce (£) } + Ng (£) (3.3) 


These three equations show that the final signal sampled 
by the Digital Programmable Oscilloscope is equal to the 
transfer function product of the amplified pulse, the 
transmitting antennas and the receiving antenna times the sum 
of the antenna coupling, the chamber response, the target 
response and its respective interaction with the chamber. For 
equation 3.2 the last two transfer functions are not 
considered since there is no target present. The total noise 


described before is added to each of these measurements. 
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C. TRANSIENT RESPONSE EVALUATION 

The transient response solution of the measured target, 
due to a double Gaussian pulse for this work, is obtained by 
an optimal deconvolution estimator developed by Raid 
(Ref. 6]. 

The first step in the process is subtracting the 
background measurement from the calibration and target 
measurements, to eliminate the antenna coupling and chamber 


clutter effects. These subtractions give the following 


equations: 

he(ft)=Y,(f) - Y;(f) (3.4) 
@ (f£)=X(f)H,(£)H,(£) (Hye(£) +H,,,(£) )} + Nz(f) (3.5) 
N,(£)=N,(£) - N5(f£) (3.6) 
meee) =i. (f£) = Y;(f) (3.7) 
ee) =X(£)HO(f)H,(f£) (HL. (£) +H. (f+ N,(f) (3.8) 
No(f)=N,(f) - N,(f) (3.9) 


The time domain subtracted measurements are transformed by 


the Fast Fourier Transform (FFT) to the frequency domain, 
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giving Y,(f) and Y,(f), and the optimal deconvolution estimator 
is applied, 
YiCf) vec 


Xo(£) Hys(f£) = —~—_—__ (Hg (£)) ys (£) 1 (3.10) 
Y,(£) Y; (one 


where X,(f)H,,(f£) 1s the target transient response due to the 
double Gaussian pulse denoted by x,(t). The smoothing 
parameter that avoids noise enhancement at frequencies where 
the spectrum of the subtracted calibration measurement 
approaches zero is denoted by "c". | 

The term X,(f)H,,(f) represents the computed calibration 
sphere response due to the specified pulse waveform 
illuminating the same size sphere as employed in the 
calibration measurement. This computed response is obtained by 
two programs written by M.A. Morgan. The first is programmed 
in Basic, called MIE, which calculates the magnitude and phase 
of the sphere's transfer function. The second, programmed in 
Fortran and called TSCT, takes the files created in the first 
and calculates the transient scattering response of the 


calibration sphere due to an incident double Gaussian pulse. 


Xy(£) = F (x(t) ) (322) 
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The double Gaussian pulse waveform used in this thesis is 


given by: 


-ay(t-tg)* @a2(t-to)” C3 .L2) 


X(t) = Cye + C2 


where C, and C, are selected to give unit value at t,: 
Xo(t>) = 1 (3.13) 


and zero d.c. content, 
[Lx (ty at = 0 (3.14) 


In equation 3.12, the double Gaussian pulse waveform is 

specified by two 10% pulse widths, At,, such that: 
oy) (3.15) 
= 0.1 for k=1,2 

For this work, the narrow pulse width was selected to be 0.15 
Nsec and the wide pulse width to be 0.30 Nsec. The reason for 
these values is that the dual parallel amplifiers, antennas 
and the DPO configuration, provide an effective 1 to 12.4 GHz 


passband. 
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IV. MEASUREMENT PROCESS AND SIGNAL PROCESSING 


A. DATA ACQUISITION 

The data acquisition in the TESL is accomplished by a 
program written in Microsoft GW Basic, called ACQU3 [Ref. 4]. 
The program begins by asking the user if it is a new 
measurement or if it is desired to use the last data. After 
this selection a friendly menu is presented on the screen 
where the user can select the filename for data output, (i.e. 
calibration, background or target). The number of time points 
in the acquired waveform is then selected. This has to be a 
power of two, ranging from 128 to 1024. The third item is the 
number of subaverages, which also has to be a power of two, 
ranging from 128 to 2048. In this thesis all the acquired 
waveforms were selected to be 1024 time points resolution and 
2048 subaverages, in a 20 nsec window, to obtain an optimal 
resolution. The fourth item is the number of data blocks in 
one record, which normally is one. The eleventh item is the 
DPO maximum vertical scale which is usually selected to be 16 
mV. The last item is the delay time of the scattering signal. 
This is a very important parameter due to the close proximity 
of the transmitting and receiving antennas. If no delay time 
were selected, the direct coupling between antennas would be 


too large, thus inhibiting scattering observance. The delay 
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time selected is usually 38 nsec. The four other items not 
mentioned help to identify the target, waveform type, date and 
hour of the measurement. When the changes have been completed, 
the user has to press the letter "R" to be able to run the 
program. After this step, the user presses the letter "B" to 
specify the bandwidth to be used. The default bandwidth is 
12.4 GHz and the program allows the user to change the default 
to 20 GHz if desired. Once the 12.4 GHz bandwidth has been 
specified, the waveform is acquired. The user has to wait 
until the DPO completes the number of subaverages selected to 
digitize the data. The acquisition and digitization for 1024 
time points and 2048 subaverages takes approximately 30 
Minutes per waveform. 

The ACQU3 program creates a file containing the raw data 
under the filename selected at the begining with extension 
WFM. Further information about the acquisition software can be 
found in Reference 4. 

Next, the user has to edit the files already created by 
the acquisition program, erase the first four lines of the 
header and rename the file with extension DAT to use the 


deconvolution algorithms developed in Matlab. 


B. DECONVOLUTION ALGORITHM 
The optimal deconvolution algorithm described in Chapter 
III section C, was initially written by M.A. Morgan for the 


Tektronis microcomputer [Ref. 2]. When the TESL configuration 
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was changed and an IBM-AT microcomputer was included, the 
algorithm was converted to Fortran 77 by Sompaee [Ref. 4]. 
In this thesis, the algorithm has been converted to Matlab 
which offers a faster and deeper analysis of the process. Two 
approaches have been developed which are described in the next 
subsections. 
1. Frequency Domain Deconvolution 

The frequency domain deconvolution algorithm (FDDA) 
uses the optimal deconvolution estimator as it was described 
in Chapter III, equation 3.10. Once it has been calculated, 
the inverse FFT is taken to find the scattered field due to a 
double Gaussian pulse. Typical waveforms of the process will 
be presented in this subsection, thus comparing the 
measurements and results obtained by two different delay lines 
for the low frequency amplifier. 

The program begins by asking the filename of the 
measured calibration sphere waveform (15 cm diameter). The 
filename entered does not required the extension DAT. Then the 
program asks for the filename of the measured background for 
the calibration sphere. After reading each of these files, the 
time window, number of time points and number of subaverages, 
Will appear on the screen. The program then plots the overlay 
of the measured waveforms, the subtracted calibration waveform 


and its spectrum as shown in Figures 7 to 10. 
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Figure 9. 
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Figure 10. Calibration Spectrum with 25.5 cm Delay Line 
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The program then asks for the filename of the measured 
target waveform and the user has the option to use the same 
background waveform used for the calibration sphere or another 
background waveform. The program presents the same information 
and plots that were presented for the calibration sphere. 
Figures 11 to 14 show the 4" thin wire 90 degrees aspect 
measurement and the target spectrum for two differents delay 
lines. After this, the program asks for the filename of the 
theoretically computed calibration sphere, which has already 
been calculated by the MIE and TSCT program as it was 
explained in Chapter III Section C, and plots the scattered 
response and its spectrum, as shown in Figure 15. 

The next step requests entry of the smoothing 
parameter "Cc". The best results in this work were found by 
using "C" = 0.1. At this point the program has already 
computed the Fast Fourier Transform (FFT) of the subtracted 
calibration waveform, the subtracted target waveform and the 
computed calibration sphere. As a part of the optimal 
deconvolution estimator, the deconvolution component, which is 
given by equation 4.1, is formed and is shown in Figures 16 


and 17 for two different delay lines. 


LOGE SieiGa) 


Z(f) = (4.1) 


Y,(£) ¥,(f) + C 
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Figure 13. 4" Thin Wire Broadside Spectrum with Minimum Delay 
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Figure 14. 4" Thin Wire Broadside Spectrum with 25.5 cm Delay 
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Theoretically Computed Calibration Sphere 
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Figure 16. 
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Figure 17. 
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The deconvolution component is then multiplied by the 
FFT of the computed calibration sphere to complete the Riad's 
optimal deconvolution estimator. The deconvolution spectrum is 
then formed, as shown in Figures 18 and 19, for the 4" thin 
wire (broadside) and for two different delay lines. After this 
step the program takes the inverse FFT and plots’ the 
deconvolved field. Figures 20 and 21 show the deconvolved 
field for the 4" thin wire (broadside) for the two different 
delay lines. The program contains a "meta" file instruction 
which saves all the plots described in the process. 

As previously explained in Chapter II Section B, the 
minimum coaxial line was set up in each amplifier. When the 
first set of measurements for the 4" thin wire were obtained 
and the data was processed by the deconvolution algorithm 
written in Fortran 77, an oscillation was noted at the front 
of the deconvolved field, as shown in Figure 20. When the 
deconvolution spectrum was compared with that obtained by 
Walsh (Ref. 5], no nulls were noted in the spectrum as were 
found in his work before he discovered the correct length for 
the delay line of the low frequency amplifier. For this 
reason, it was thought that there was not a phasing problem 
between the low and high frequency amplifiers. Several 
measurements were done for canonical targets and all the 


deconvolved fields showed the same undesired oscillation. 
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The first attempt at elimination of the oscillation 
was directed towards reduction of the quantization of the DPO. 
Four measurements were required to do this. A calibration 
measurement with its respective background was performed using 
the DPO maximum vertical scale set to 16 mV and the 4" thin 
wire (broadside) with its background set to 8 mV maximum 
vertical scale. The same procedure was followed for other 
canonical targets like the 8 cm and 12 cm diameter spheres 
without improvement in the deconvolved fields. 

The second trial was the modification of the 
deconvolution algorithm written in Fortran 77, to save the 
subtracted calibration and target waveforms. These files were 
then edited and the subtracted waveforms modified. The points 
before and after the response were zero padded. When the 
modified data were processed the deconvolved field still had 
the oscillation. The deconvolution algorithm was converted to 
Matlab and the same previous trials were performed with nearly 
identical results to those found using Fortran. A third trial 
modified the deconvolution component shown in Figures 16 and 
17. The d.c. component, the frequency components up to 1 GHz 
and the frequency components from approximately 13 GHz to 25 
GHz were made zero and it was multiplied by the FFT of the 
computed calibration sphere, to obtain the deconvolution 
spectrum. When the inverse FFT was taken to obtain the 
deconvolved field, no improvement was noted. The fourth trial 


included the measurements taken with an extra piece of a low 
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density styrofoam pedestal shown in Figure 1, because it 
looked like the targets were not perfectly aligned with the 
center of the antennas panel. Again, no improvement was 
obtained when the measurements were processed. 

Next, the deconvolution algorithm written in Matlab 
was modified and called Time Domain Deconvolution. This is 
explained in more detail in the next subsection. Equation 4.1, 
which represents the deconvolution component in the frequency 
domain, was transformed to the time domain and the deconvolved 
impulse response was formed. This impulse response was 
convolved with the time domain computed calibration sphere to 
form the hybrid deconvolved field. The same undesired response 
was observed. 

During another test the deconvolved impulse response 
waveform shown in Figure 24 was modified. The points before 
the impulse response were made zero. When this modified 
waveform was convolved with the time domain computed 
calibration sphere, the oscillation at the front of the hybrid 
deconvolved field disappeared. It was thought that the 
oscillation could be caused by a d.c. component present in the 
measurements. A new set of complete measurements were done 
with a d.c. blocking capacitor installed at the output of the 
receiving antenna. After the signal processing the oscillation 
still existed. 

A new program in Matlab was created to verify the 


effect of the delay line in both amplifiers. A set of three 
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measurements required by the process were obtained with just 
the low frequency amplifier while the other amplifier was 
replaced by a dummy load. The same procedure was followed for 
the high frequency amplifier. The program begins by asking the 
user how many time points to shift the low frequency signal. 
The program then asks for the filenames of each of the six 
waveforms and forms four subtracted waveforms as it was 
explained for the deconvolution algorithm. Once this has been 
done the program adds the corresponding subtracted waveforms 
for both amplifiers and follows the procedure described 
before. It was determined, in case of the values shown in 
Chapter II Section B, that each time point is equivalent to 
approximately 4 mm of delay line. Several trials were 
performed, sweeping from 1 to 100 time points for the low 
frequency Signals. It was noted that the delay line had a 
large effect in the final deconvolved field. The same 
procedure was performed to shift the high frequency signals. 
The user must keep in mind that the program has been designed 
to shift just the low frequency signals, so in order to shift 
the high frequency signals, the user must input the data 
filename of the 6-12 GHz when the program requests the data 
filenames for the 1-6 GHz and vice versa. 

The best result was obtained by shifting the signal of 
the low frequency amplifier by 30 time points. These 30 time 
points are equivalent to an extra delay of 12 cm, as shown in 


Chapter II Section B. These measurements were done with a 
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minimum delay line for both amplifiers. With this result, a 
25.5 cm delay line was constructed for the low frequency 
amplifier, thus giving 13.5 cm corresponding to the minimum 
delay and 12 cm corresponding to the extra delay. When the 
measurements were done with the new delay the undesired 
oscillation was not present, as is shown in Figures 21 and 27. 
Several sets of measurements for canonical targets were done 
with this new delay line yielding excellent results that agree 
with the theoretical waveforms, as will be shown in the next 
section. 

It can be clearly observed in Figures 7 and 11, in 
comparison with Figures 8 and 12 presented in this subsection, 
that with the minimum delay line the low frequency bandpass 
Signal was arriving at the DPO sampling head before the high 
frequency component. After the 25.5 cm delay line was 
installed, both signals are arriving at the DPO sampling head 
at the same time. The effect of the delay line in the final 
deconvolved field is critical to obtain desired waveforms, as 
shown in Figure 21. Conversely, if the delay line does not 
provide simultaneous signal arrival, an unacceptable waveforn, 
which includes a series of additional oscillations, will 
occur. A simulation using a 35.5 cm delay line (55 time 
points) was introduced for the low frequency amplifier 
yielding Figure 22. This is a very undesirable waveform 
although its spectrum, Figure 23, does not display any deep 


nulls as presented by Walsh (Ref. 5]. 
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Figure 22. 


Simulated 35.5 cm Delay Line 
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Figure 23. 4" Thin Wire Broadside Deconvolution Spectrum with 
a Simulated 35.5 cm Delay Line 


2. Time Domain Deconvolution 

The time domain deconvolution algorithm (TDDA) is a 
continuation of the frequency domain deconvolution algorithm 
(FDDA). The optimal deconvolution estimator described in 
Chapter III can be evaluated completely in the time domain for 
the case of no smoothing (e.g. C=0). 

Both the FDDA and TDDA were computed in the same 
program. The FDDA calculates the deconvolved field by 
multipliying the deconvolution component, presented in 
equation 4.1, with the FFT of the computed calibration sphere 
followed by an inverse FFT. Unlike the FDDA, the TDDA 
calculates the deconvolved field by taking the inverse FFT of 
the deconvolution component, forming the deconvolved impulse 
response, and convolving it with the computed calibration 
sphere as shown in equation 4.2. 

Mott) hye ( yer oaks *& [X_(t) hys(t) J (4.2) 
Wert) Y,(E) + € 

The two components of the convolution specified by 
equation 4.2, the computed calibration sphere and the 
deconvolved impulse response, are shown in Figures 24 and 25 
for a 4" thin wire (broadside) and for two differents delay 
lines. The hybrid deconvolved fields are shown in Figures 26 
and 27. Figure 27 displays excellent agreement with its 


respective theoretical one. 
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Figure 25. Computed Calibration Sphere and 4" Thin Wire 


Broadside Deconvolved Impulse Response with 25.5 cm Delay Line 
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Figure 26. 4" Thin Wire Broadside Hybrid Deconvolved Field 
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with 25.5 cm Delay Line 


50 


C. SYSTEM VALIDATION 

Validation of the upgraded Transient Electromagnetic 
Scattering Laboratory was performed by comparing 
experimentally obtained deconvolved fields with theoretically 
computed waveforms for canonical targets. These canonical 
targets are: 

1) 8 centimeter diameter metal sphere. 

2) 12 centimeter diameter metal sphere. 

3) 10 cm long, 2.36 millimeter diameter thin wire 90 

degrees aspect (broadside). 

4) 10 cm long, 2.36 millimeter diameter thin wire 30 

degrees aspect. 

The theoretical responses for the 8 cm and 12 cm spheres 
are obtained by the two programs, MIE and TSCT, as mentioned 
in Chapter III. MIE calculates the magnitude and phase of the 
sphere's transfer function. The inputs used in this program 
are summarized in Table 2. TSCT takes the files created by MIE 
and computes the scattering response due to an incident double 
Gaussian pulse. Table 3 summarizes the inputs used in the TSCT 
program. 

Another program written by M.A. Morgan, called TDIE-DG, 
computes the theoretical scattering response for the thin wire 
due to a double Gaussian pulse. The program makes the 
discretization and solves a time domain integral equation to 
calculate the induced currents in time and space. The program 


then computes the backscattering far field, using the 


ol 


currents, by a numerical time space integration. The user 


supplied inputs used in TDIE-DG are summarized in Table 4. 


TABLE 2. NUMERICAL COMPUTATIONS INPUTS OF MIE PROGRAM 









————————— eS EE oe 


Enter filename of magnitude and Selected by the 

















phase array user 


Sphere radius in meters 0.04, 0.06 
oo <a 





a 


TABLE 3. NUMERICAL COMPUTATIONS INPUTS OF TSCT PROGRAM 



















As selected by the 








Enter data filename for magnitude 









and phase created by MIE user in MIE 


| Narrow 10% PW in nsec Om 5 
Wide 10% PW in nsec 0.30 
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TABLE 4. NUMERICAL COMPUTATIONS INPUTS OF TDIE-DG PROGRAM 


a 
a 
LS Oe oO 












0.00117 








Far field computation 
pene ey 


| Theta in degrees for far field 





| computation = 180 - Alpha 





Figures 28, 29, 30, and 31 compare the theoretical and 
deconvolved waveforms for the canonical targets described in 
this section. The high level of accuracy obtained by the new 
TESL configuration and the algorithms used are apparent from 


these examples. 
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Figure 28. 
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Figure 30. 
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D. NOISE PERFORMANCE AND SPECTRUM 

The noise performance of the TESL was evaluated by 
comparing the root mean square (RMS) value of the noise 
waveforms from old and new configurations. Noise waveforms are 
formed by subtracting two consecutive background measurements. 
Although there is record of the old configuration background 
measurements by Walsh [{Ref. 5], two consecutive measurements 
were not available. Because of this, the noise performance of 
the updated TESL was calculated using only one background 
measurement subtracted from the calibration sphere measurement 
so that the old and new configuration noise performance could 
be compared during the final 10 Nsec or so of the sampling 
window, where only noise exists. 

Figure 32 shows both subtracted calibration waveforms. 
This particular window is enlarged and shown in detail in 
Figure 33. 

A Matlab program was written to perform this procedure and 
calculate the RMS value of both noise waveforms. Due to the 
limited available data for the old configuration, only a few 
comparisons could be done. In all cases the new TESL 
configuration resulted in better noise performance. Figure 33 
shows the new configuration RMS value of 0.024 mV compared to 


the previous setup which yields 0.034 mV. 
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Figure 33. 


Old and New Configuration Noise Waveforms 
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Another Matlab program was written to calculate and plot 
the noise spectrum. The program just requires two consecutive 
background measurements and asks for the smoothing factor, 
which is entered by the operator. When a factor of 1 was 
utilized, the unsmoothed noise spectrum in Figure 34 was 
obtained. Figure 35 shows the time-smoothed system noise 


spectrum when a factor of 8 is used. 
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Figure 34. Unsmoothed Noise Spectrum 
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Figure 35. Time Smoothed Noise Spectrum 
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E. SIGNAL TO NOISE RATIO 

A Matlab program was written to calculate the SNR of the 
TESL. The program begins by asking for the filenames of the 
measured target waveform and its background. The signal plus 
noise (denoted by s(t,)) is formed by subtracting this 
background from the target waveforn. 

Next, the program asks for the two consecutive background 
waveforms to estimate the system noise. The noise waveform 
(denoted by n(t,)) is then formed by subtracting one of these 
background waveforms from the another. After this, the program 
presents a message on the screen asking the user to view the 
record to select the desired time-window which contains only 
the background signal. When the beginning and the end of the 
time window have been selected, the program calculates the SNR 


as shown in equation 4.3. 


y (s*(t,) - n°(e)] 


SNR = 10 Log [ OTe 
k 


} (4.3) 


where s(t,) 1s the signal plus noise waveforn, 
n(t,) is the noise waveform. 

McDaniel, Sompaee, and Walsh (Refs. 3,4 & 5] calculated 
the SNR using a 10 cm long and 5.2 mm diameter thick wire as 
a target in the same way as desribed in this section. Walsh 
[Ref. 5] reported a SNR of 26.6 dB with this target. This 
target was not available when this thesis started. A 10 cm 


long and 2.36 mm diameter thin wire was used as a target in 


64 


order to estimate a comparison with the SNR obtained by the 
previous work. The same number of points, same number of 
subaverages and the same acquisition software were used. 
Although this target is thinner than that used before and a 
smaller response was expected, the SNR obtained was 30.56 GB 
using a time-window from 1.9 nsec to 7 nsec, as shown in 
Figure 36. Different sets of measurements of this thin wire 
were analyzed, yielding approximately the same SNR. These 
calculations show that a significant improvement has been 


achieved in the SNR of the TESL by the results of this work. 


F. TARGET LIBRARY 

The high level of accuracy of the upgraded TESL has been 
shown in Section C of this Chapter for simple canonical 
targets through comparisons with computations. These canonical 
target validations indicate the accuracy to be expected for 
complex targets, such as tactical aircraft models, where 
computations are exceedingly difficult. 

A target library was created to support research in aspect 
invariant radar target identification based on natural 
resonances. This library contains six files, each one 
corresponding to a different 1/72 scale model military 
aircraft which has been coated with silver paint. Each file 
contains 9 measurements. Each measurement was acquired with 


the same parameters described in Section A of this Chapter. 
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Figure 36. Signal and Noise Waveforms 
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(AW) painscayy 


The 9 measurements for each aircraft are: 
- A calibration sphere measurement 
first background measurement 
target measurement (0 degrees) 
second background measurement 
target measurement (30 degrees) 
third background measurement 


target measurement (90 degrees) 


: 
rF FP FP FP FP YP PY 


fourth background measurement 
- A target measurement (180 degrees) 
The full size dimensions of the aircraft selected — 
specified in Table 5. Figures 40 and 41 show some of the 
waveforms obtained. 


TABLE 5. AIRCRAFT FULL SIZE DIMENSIONS 


Length 19.20 | 16.82 18.85 | 20.84 

(meters) 

Height 3502 4.32 
(meters) 

Wingspan | 13.10 | 14.26 | 10.80 | 11.52 10.60 
(meters) 

- i ee 


67 












Di<15 
F] 
ia! 












Tailplane 





(meters) 





See = | H 


-I\~r i\ pel 





(sasu) auuty, : 


eae | 
Cay 





S| ee 


Pov JeIV UORV 
(oasu) awe 
I e iG I 
i 


Ai ee, 


TAVUOULV 











(w/Au) (1)F 


(w/Aw) (1a 


(aasu) atany 


C l 


¢ b ¢ 
. ss : 
| 
| 
= | 
| | 
| | | 
| ! ; 
| 
: 
ee Spe. * 7 
€C # JAVUOUV 
(oasu) aut, 
g b t v4 
: | 
pa 
| | 
: | 3 
| | ; 
| | 
| : 
i ee ee 


L # LAV UOUIV 


0) 


i——-}- 


-|Z- 





(w/Am) ()q 


(my/Atm) (3)g 
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V. CONCLUSIONS 


A. SUMMARY 

This thesis describes several improvements obtained by 
the upgraded TESL. The first stage of this effort was the 
replacement of the old HP 8349A 2-20 GHz microwave 
preamplifier with a new Avantek AWT 13533 5-13 GHz (which has 
a higher gain and a lower noise figure than the HP amplifier). 
A minimum length coaxial line was used with each amplifier. An 
undesired oscillation in front of the deconvolved scattered 
Signal for different measured targets was found. The first 
approach to eliminate the oscillation was directed towards 
reduction of the quantization of the DPO, without improvement 
in the final response. Another test was performed taking a set 
of measurements with a d.c. blocking capacitor installed at 
the output of the receiving antenna. After deconvolution 
Signal processing, the same undesired response was observed. 

A Matlab program was written to evaluate the effect of the 
delay line in each amplifier. It was found that the delay 
lines had a large effect in the deconvolved signal. The best 
results were obtained when the measured waveform from the 2-6 
GHz amplifier was delayed 30 time points. It was determined 
that each shifted time point is equivalent to 4 mm, so an 


extra 12 cm length delay line for the low frequency amplifier 


70 


was needed. With this result, a 25.5 cm delay line was 
constructed for the 2-6 GHz amplifier, thus giving 13.5 cm 
corresponding to the minimum delay and 12 cm corresponding to 
the extra delay. After sets of measurements were taken and 
processed, the deconvolved signals were validated by comparing 
with those theoretically predicted. The updated TESL was shown 
to yield a high level of accuracy. 

Additionally, Matlab programs to calculate the SNR, noise 
spectrum and performance were written and tested. An 
improvement in the SNR of 4 a@B was found, although the 
reference target was similar but of thinner diameter than the 
one used in previous experiments. 

The final stage in this thesis was the compilation of a 
high fidelity library for canonical and complex targets. This 
library will support future research into aspect invariant 
radar target identification and other scattering 


experimentation. 


B. FUTURE CONSIDERATIONS 
Increased SNR can be obtained by replacing the one-watt 
amplifiers (1-6 GHz and 6-12 GHz) with two-watt amplifiers. 
Another improvement to the TESL configuration would be 
mounting all of the hardware in a vertical system rack which 
will allow higher system flexibility for future modifications. 
Since there is a time domain radar target library then it 


would be worthwhile to create a corresponding frequency domain 


Taek 


library to facilitate complete TESL experimental analysis. 
Presently there is an ongoing effort in this area which 


employs a stepped-frequency continuous wave system based on an 


HP 8510 network analyzer. 
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APPENDIX A. TIME AND FREQUENCY DOMAIN DECONVOLUTION PROGRAM 
The following Matlab program calculates the Deconvolved 

Field using the Raid's optimal deconvolution estimator as 

described in chapter IV. 

NAME OF THE PROGRAM: DECON.M 

ALGORITHM MADE BY Dr. M.A. MORGAN. 

WRITTEN BY Lt. ALDO E. BRESANI. 


DECONVOLUTION ALGORITM FOR THE NPS TRANSIENT ELECROMAGNETIC 
SCATTERING LABORATORY. 


INPUTS ARE THE MEASURED WAVEFORMS, COMPUTED CALIBRATION 
SPHERE AND THE SMOOTHING PARAMETER. 


RFF FF KKE 


clg 
Clear 
filename=input(’Enter the ’’Calibration filename’’.dat to load :’,’s’); 
eval([’load ’,filename,’.dat’)) 
datal=eval(filename); 
nsavl=datal(1); % NUMBER OF SUBAVERAGES 
ntpl =datal(4); % NUMBER OF TIME POINTS 
tstart1=0; 
tend 1=datal(2); 
tl=tend1-tstartl; 
[Time window is ’,num2str(tl ),”’ nanosec’] 
[Number of Time Points is ’,num2str(ntp1),’’] 
[Number of subaverages is ’,snum2str(nsavl ),’’] 
cal=zeros(I ,ntp1); 
for j=l:ntpl 

cal(j)=datal (j+4); 
end ; 


filename=input(’Enter the ’Calibr.Background filename’’.dat to load :’,’s’); 
eval([’load ’,filename,’.dat’}) 
data2=eval(filename); 

nsav2=data2(1); % Number of Subaverages 
ntp2=data2(4); % Number of Time Points 
tstart2=0; 

tend2=data2(2); 

t2=tend2-tstart2; 

[Time window is ’,numd2str(t2),’ nanosec’] 
("Number of Time Points is ’,num2str(ntp2),’’] 
("Number of subaverages is ’,num2str(nsav2),’’} 
back 1=zeros(1,ntp2); 
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for j=I:ntp2 
back 1(j)=data2(j+4); 
end ; 
ntp=ntp]l; 
dt=tl/(ntp-1); 
t=0:dt:tend]; % DEFINING TIME AXIS 
f=((ntp/(2*t1))/ntp)*(O:ntp-1); % DEFINING FREQUENCY AXIS 
dt2=10/511; 
(2=0:0t2 10: 
dtl=5/255; 
T1=0:dt1:5; 
ycal=zeros(1,ntp); 
for k=I:ntp 
ycal(k)=cal(k)-back 1(k); 
end; 


subplot(211),plot(t2,(1000.*(cal(1:512))),’-r’,...... 

‘agen t2,(1000.*(back1(1:512))),’--g’) 

title’ CALIBRATION AND BACKGROUND WAVEFORM’); 
xlabel(’Time (nsec)’); 

grid; 

ylabel(’Measured (mV)’); 

pause; 


subplot(212),plot(t2,(1000.*(ycal(1:512)))) 
title’/SUBTRACTED CALIBRATION WAVEFORM’); 
xlabel(’Time (nsec)”’); 

grid; 

ylabel(*7Measured (mV)’); 

meta lica; 

pause; 

clg; 


Yeal = fft(ycal,2*ntp); 


plot(f,(1000*%abs(Ycal(1:ntp)))) 

title”? CALIBRATION SPECTRUM’); 
xlabel( Frequency (GHz)’); 

grid; 

ylabel(’Magnitude’); 

meta lica; 

pause; 

clg ; 


filename=input(’Enter the Target filename’’.dat to load :’,’s’); 
eval([load ’,filename,’.dat’]) 

data3=eval(filename); 

nsav3=data3(1); % Number of Subaverages 

ntp3=data3(4); % Number of Time Points 

tstart3=0; 

tend3=data3(2); 
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t3=tend3-tstart3; 
("Time window is ’,num2str(t3),’ nanosec’} 
("Number of Time Points is ’,num2str(ntp3),’} 
[’Number of subaverages is ’,snum2str(nsav3),”"} 
targ=zeros(1 ,ntp3); 
fer j=i-ntp3 

targ(j)=data3(j+4); 
end ; 


h=input(’Do you want to use the same Back. file (0 for N, | for Y):’), 
if h== 


filename=input(’Enter the ’ Targ.Background filename”’.dat to load.... ....:’,’s’); 
eval([(’load ’,filename,’.dat’}) 
dataS5=eval(filename); 
nsav5=data5(1); % Number of Subaverages 
ntpS=data5(4); % Number of Time Points 
tstart5=0; 
tend5=data5(2); 
t5=tend5-tstart5; 
(Time window is ’,num2str(t5),’ nanosec’} 
{Number of Time Points is ’,num2str(ntp5),”’} 
{Number of subaverages is ’,num2str(nsav5),’’]} 
back2=zeros(1,ntp5); 
for j=I:ntp5 

back2())=data5(j+4); 
end ; 


else 
back2=zeros(1,ntp); 
for j=I:ntp2 
back2(j)=back1(j); 
end; 
end; 


ysct=zeros(1:ntp); 

for k=I:ntp 
ysct(k)=targ(k)-back2(k); 

end; 


subplot(211),plot(t2,(1000.*(targ(1:512))),’-r’,...... 

a t2,(1000.*(back2(1:512))),’--g’) 

title? TARGET AND BACKGROUND WAVEFORM’), 
xlabel(’Time (nsec)’); 

ylabel( Measured (mV)’); 

grid; 

pause; 


subplot(212),plot(t2,(1000.*(ysct(1:512)))) 
title’SUBTRACTED TARGET WAVEFORM’); 


Ue 


xlabel(’Time (nsec)’); 
ylabel(’ Measured (mV)’); 
grid; 

meta lica; 

pause; 

Cie 


Ysct = fft(ysct,2*ntp); 


plot(f,(1000*abs( Ysct(1:ntp)))) 
title’ Target Spectrum’): 
xlabel(’ Frequency (GHz)’); 
ylabel(’Magnitude’); 

grid; 

meta lica; 

pause; 

cle ; 


filename=input(’Enter the ’Computed Sphere filename’’.dat to load.... .... Ss 
eval((load ’,filename,’.dat’]) 
data4=eval(filename); 
ntp4=data4(1); % Number of Time Points 
tstart4=data4(2): 
tend4=data4(3): 
t4=tend4-tstart4; 
(Time window is ’,num2str(t4),’ nanosec’] 
("Number of Time Points is ’,num2str(ntp4),’’] 
ycom=zeros(! ,ntp4); 
for j=l:ntp4 

ycom(})=data4(j+3); 
end ; 


l=input(’Enter the smoothing factor :’); 


YIcal= conj(Ycal) ; 
a=0; 
for n=1:2*ntp 
a=a+(abs( Ycal(n))).*2 ; 
end; 
b=a/(ntp); 


for n=1:(2*ntp) 

Z(n)=(Ysct(n).* YI cal(n))./((abs( Ycal(n)))*2+d);:% DECOVOLUTION COMPONENT 
end; 

z=zeros(1,2*ntp) ; 

z=ifft(Z,2*ntp) ; % DECONVOLVED IMPULSE RESPONSE 
Ycom=fft(ycom,2*ntp) ; 
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s3=Ycom.*conj(Ycom) ; % SPECTRUM COMPUTED CALIBRATION SPHERE 
subplot(211),plot(T1,(1000.*ycom(1:256))) ; grid 

title(*COMPUTED SPHERE SCATTERED E(t) IN mV/m’); 

xlabel(’Time (nsec)”’); 

pause; 


subplot(212), plot(f,(1000.*real(s3(1:ntp)))) 

grid: 

title COMPUTED SPHERE SCATTERED SPECTRUM |E(f))’): 
xlabel(’ Frequency (GHz)’); 

meta lica; 


Ziel.*Ycom - 
% FREQUENCY DOMAIN DECONVOLUTION 
zz=ifft(ZZ,2*ntp); 


plot(f,(1000.*abs(Z(I:ntp)))) ; 

title’? DECONVOLUTION COMPONENT?); 
grid; 

xlabel( Frequency (GHz)’); 

meta lica; 

pause; 

clg 


plot(f,(1000.*abs(ZZ(1:ntp)))) 

title’ DECONVOLUTION SPECTRUM, SMOOTHING: ’,num2str(1)J); 
xlabel( Frequency (GHz)’): 

ylabel(’ Magnitude’); 

grid; 

meta lica; 

pause; 

clg 


plot(T1,(1000.*real(zz(1:256)))) 

title(/ DECONVOLVED FIELD, SMOOTHING: ’,num2str(1)}); 
xlabel(’Time (nsec)’); 

ylabel(’E(t) (mV/m)’); 

grid; 

meta lica; 

pause; 

cle 


% TIME DOMAIN DECONVOLUTION 
cc=conv(z(1:256),ycom(1:256)); 


subplot(211),plot(T1 ,(1000.*ycom(1:256))) ; grid 
title’ COMPUTED SPHERE SCATTERED E(t) IN mV/m’); 


val 


xlabel(’Time (nsec)’); 

pause; 

subplot(212),plot(T1 ,(1000.*real(z(1:256)))) ; 

title(’7 DECONVOLVED IMPULSE RESPONSE, SMOOTHING: ’,numé2str(1))); 
xlabel(’ Time (nsec)’); 

grid; 

meta lica; 

pause; 

clg; 


plot(T1,(1000.*(cce(1:256)))) ; 

title(’ HYBRID DECONVOLVED FIELD, SMOOTHING: ’,num2str(1)}); 
xlabel(’ TIME (nsec)’); 

ylabel(SCAT.E-FIELD(mV/m)’); 

grid; 

meta lica; 

pause; 

end; 
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APPENDIX B. DELAY LINE SIMULATION PROGRAM 
The following Matlab program simulates the effect of the 


delay line in each of the amplifiers of the TESL. 


% NAME OF THE PROGRAM: SIMDEL.M 
% WRITTEN BY Lt. ALDO E. BRESANI 


% INPUTS ARE THE MEASURED WAVEFORMS, COMPUTED CALIBRATION 
SPHERE AND THE SMOOTHING PARAMETER. 


cle 
clear 


q=input(’Enter the shift for 1-6 GHz: ’); 


filename=input(’Enter the ’Calibration 1-6 GHz’’.dat to load :’,’s’); 
eval([’load ’,filename,’.dat’]) 
datal=eval(filename); 
ntp! =datal(4); % NUMBER OF TIME POINTS 
tstart!=0; 
tend! =data!(2): 
tl=tend]-tstart]; 
call=zeros(1,ntp!); 
for j=I:ntp! 
call (j)=datal(j+4); 
end ; 
call ]=zeros(1:1024); 
for k=1:q 
call 1(k)=0; 
end; 
for k=(q+1):1024 
call 1(k)= call(k-q); 
end; 


filename=input(’Enter the *’Calibr.Background !1-6 GHz’’.dat to load :’,’s’); 
eval([’load ’,filename,’.dat’}) 
data2=eval(filename); 
ntp2=data2(4); % Number of Time Points 
tstart2=0; 
tend2=data2(2):; 
t2=tend2-tstart2; 
back ]=zeros(1 ,ntp2); 
for j=I:ntp2 
back 1(j)=data2(j+4); © 
end ; 


Ue 


back 1 1=zeros(1:1024); 
for k=l:q 
back11(k)=0; 
end; 
for k=q+1:1024 
back1 1(k)= back1(k-q); 
end; 


ntp=ntpl: 

dt=tl/(ntp-1); 

t=0:dt:tend 1; % DEFINING TIME AXIS 
f=((ntp/(tl))/ntp)*(O:ntp/2-1); % DEFINING FREQUENCY AXIS 
dt2=10/511; 

t2=0:dt2:10; 

dtl=5/255: 

T1=0:dtl:5; 


ycall=zeros(1,ntp); 
for k=I:ntp 

ycall(k)=call 1(k)-back1 1(k); 
end; 


subplot(211),plot(t2,(1000.*(cal11(1:512))),’-r’,...... 

sae t2,(1000.*(back11(1:512))),’--g’) 

title °CALIBRATION AND BACKGROUND WAVEFORM’); 
xlabel(’Time (nsec)”’); 

ylabel(’Measured (mV)”); 

grid; 

pause; 


subplot(212),plot(t2,(1000.*(ycall(1:512)))) 

title (SUBTRACTED CALIBRATION WAVEFORM’); 
xlabel(’Time (nsec)’); 

ylabel( Measured (mV)’): 

grid; 

pause; 

clg; 


filename=input(’Enter the ’’Calibration 6-12 GHz’’.dat to load :’,’s’); 
eval([load ’,filename,’.dat’]) 
data10=eval(filename); 
ntp10 =datal0(4); % NUMBER OF TIME POINTS 
tstart10=0; 
tend10=datal0(2); 
tlO=tend10-tstart10; 
cal2=zeros(1,ntp10); 
for j=I:ntp10 
cal2(j)=data10(j+4); 
end ; 


filename=input(’Enter the ’Calibr.Background 6-12 GHz’’.dat to load... s 
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eval([’load ’,filename,’.dat’]) 
data20=eval(filename); 
ntp20=data20(4); % Number of Time Points 
tstart20=0: 
tend20=data20(2); 
t20=tend20-tstart20; 
back2=zeros(1,ntp20); 
for j=I:ntp20 
back2(j)=data20(j+4); 
end ; 


ycall2=zeros(1:ntp); 

for k=I:ntp 
ycall2(k)=cal2(k)-back2(k); 

end; 


subplot(211),plot(t2,(1000.*(cal2(1:512))),’-r’,...... 

oe t2,(1000.*(back2(1:512))),’--g’) 

title” CALIBRATION AND BACKGROUND 6-12’); 
xlabel(’Time (nsec)’); 

ylabel( Measured (mV)’): 

grid; 

pause; 


subplot(212),plot(t2,(1000.*( ycall2(1:512)))) 

title (SUBTRACTED CALIBRATION WAVEFORM’); 
xlabel(’Time (nsec)’); 

ylabel(’Measured (mV)’); 

grid: 

pause; 

clg; 


for k=1:1024 
ycal(k)=ycall(k)+ycall2(k); 
end; 
for n=513:1024 
ycal(n)=0; 
end; 
Ycal = fft(ycal,ntp); 


plot(f,(1000*abs(Ycal(1:ntp/2)))) 
title(’Calibration Spectrum’); 
xlabel(’Frequency (GHz)’); 

grid; 

pause; 

clg ; 


filename=input(’Enter the Target 1-6 GHz’’.dat to load :’,’s’); 
eval({‘load ’,filename,’.dat’}) 

data3=eval(filename); 

ntp3=data3(4); % Number of Time Points 
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tstart3=0; 
tend3=data3(2); 
t3=tend3-tstart3; 
targ]=zeros(1,ntp3); 
for j=I:ntp3 

targ1(j)=data3(j+4); 
end ; 


targ]]=zeros(1:1024); 


for k=]:q 
targ] 1(k)=0; 
end; 


for k=q+1:1024 
targ11(k)= targ1(k-q); 
end; 


filename=input(’Enter the Target 6-12 GHz’’.dat to load :’,’s’); 
eval([’load ’,filename,’.dat’]) 
data30=eval(filename)- 
ntp30=data30(4); % Number of Time Points 
tstart30=0; 
tend30=data30(2); 
t30=tend30-tstart30; 
targ2=zeros(1,ntp30); 
for j=]:ntp30 
targ2(j)=data30(j+4): 
end ; 


ysctl=zeros(1:ntp); 
ysct2=zeros(1:ntp); 


for k=I:ntp 
ysctl(k)=targ11(k)-back]1 1(k); 
ysct2(k)=targ2(k)- back2(k); 
end; 


ysct=zeros( I:ntp); 


for k=I:ntp 
ysct(k)=ysct1(k)+ysct2(k); 
end; 


subplot(211),plot(t2,(1000.*(targ 1 1(1:512))),’-r’,...... 
ae t2,(1000.*(back11(1:512))),’--g’) 
title?TARGET AND BACKGROUND 1-6 GHz’); 
xlabel(’Time (nsec)’); 

ylabel(’ Measured (mY)’): 

grid; 

pause; 
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subplot(212), plot(t2,(1000.*(ysct1(1:512)))) 
title’ SUBTRACTED TARGET WAVEFORM’): 
xlabel(’Time (nsec)’); 

ylabel(’Measured (mV)’):; 

grid; 

pause; 

cle; 


subplot(211),plot(t2,(1000.*(targ2(1:512))),’-r’,...... 
ess. t2,(1000.*(back2(1:512))),’--g’) 

title( TARGET AND BACKGROUND 6-12 GHz’): 
xlabel(’Time (nsec)’); 

ylabel(’Measured (mV)’); 

grid; 

pause; 


subplot(212),plot(t2,(1000.*(ysct2(1:512)))) 

title’ SUBTRACTED TARGET WAVEFORM’); 
xlabel(’Time (nsec)’): 

ylabel(’Measured (mV)’): 

grid; 

pause; 

clg; 


Ysct = fft(ysct,ntp); 


plot(f,(1000*abs( Ysct(1:ntp/2)))) 
title(’Target Spectrum’); 
xlabel(’Frequency (GHz)”’); 

grid; 

pause; 

clg ; 


filename=input(’Enter the ’Computed Sphere filename’’.dat to load... ....:’,’s’); 
eval([{’load ’,filename,’.dat’)) 
data4=eval(filename); 
ntp4=data4(1); % Number of Time Points 
tstart4=data4(2); 
tend4=data4(3); 
t4=tend4-tstart4; 
(Time window is ’,num2str(t4),’ nanosec’] 
("Number of Time Points is ’,num2str(ntp4),’’] 
ycom=zeros(1,ntp4); 
for j=I:ntp4 

ycom(j)=data4(j+3); 
end ; 


cc=conv(ysct,ycom); 

gg=zeros(1,ntp); 

for n=1:1023 
ge(n)=cc(n); 


83 


end; 
g2(1024)=0; 
GG=fft(gg,ntp), 


l=input(Enter the smoothing factor :’); 
Ylcal= conj(Yeal) ; 
a=0 
for f=l-nte 
a=a+(abs( Ycal(n))).*2 ; 
end; 
b=a/(ntp); 


C= Z 
d=cab: 
Z=zeros(1,ntp) ; 
for n=I:ntp 
Z(n)=(GG(n).* Ylcal(n))./((abs(Ycal(n)))*2+d) ; 
end; 


plot(f,(1000.*abs(Z(1:ntp/2)))) 


title(’?7 DECONVOLUTION SPECTRUM, SMOOTHING: ’,num2str(1)}) 


xlabel(’ Frequency (GHz)’); 
ylabel(’Magnitude’); 

grid; 

meta sim; 

pause; 

clg 


z=zeros(1,ntp) ; 
z=ifft(Z,ntp) ; 


plot(T 1,real(z(1:256))) 


title(? DECONVOLVED FIELD, SMOOTHING: ’,num2str(1)}); 


xlabel(’ Time (nsec)’); 
ylabel( E(t) (mV/m)’); 
grid; 

meta sim; 

pause; 
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APPENDIX C. SNR AND NOISE SPECTRUM PROGRAM 

The following Matlab program calculates the signal to 
noise ratio and the noise spectrum of the TESL. 
NAME OF THE PROGRAM: SNR.M 
WRITTEN BY Lt. ALDO E. BRESANI 
THIS PROGRAM CALCULATES THE SIGNAL TO NOISE RATIO FOR THE 
TRANSIENT ELECTROMAGNETIC SCATTERING LABORATORY, AND THE 
SYSTEM NOISE SPECTRUM. 


INPUTS ARE THE TARGET MEASUREMENT AND ITS BACKGROUND, AND 
TWO CONSECUTIVE BACKGROUND MEASUREMENTS. 


FF FFF KK 


clg 
clear 
filename=input(’Enter the ’’Target filename”.dat to load :’,’s’); 
eval([’load ’,filename,’.dat’]) 
datal=eval(filename); 
nsavl=datal(1); % NUMBER OF SUBAVERAGES 
ntpl =datal(4); % NUMBER OF TIME POINTS 
tstart1=0; 
tend1=datal (2); 
tl=tend1-tstartl; 
(’Time window is ’,num2str(tl),’ nanosec’] 
("Number of Time Points is ’,num2str(ntp!1),”’] 
("Number of subaverages is ’,num2str(nsav1),’’] 
targ=zeros(1,ntp1 ); 
for j=I:ntp1 

targ())=data1(j+4); 
end; 


filename=input(’Enter the ’’ Target Backg. filename”’.dat to load :’,’s’); 
eval({’load ’ ,filename,’.dat’)) 
data4=eval(filename); 

nsav4=data4(1); % Number of Subaverages 
ntp4=data4(4); | % Number of Time Points 
tstart4=0; 

tend4=data4(2); 

t4=tend4-tstart4; 

(Time window is ’,num2str(t4),’ nanosec’] 
[’Number of Time Points is ’,snum2str(ntp4),”] 
["Number of subaverages is ’,snum2str(nsav4),”’] 
back0=zeros(1,ntp4); 

for j=l:ntp4 
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back0(j)=data4(j+4); 
end ; 


filename=input(’Enter the Background! filename” .dat to load :;’,’s’); 
eval({load °,filename,’.dat’]) 
data2=eval(filename); 
nsav2=data2(l1); % Number of Subaverages 
ntp2=data2(4); % Number of Time Points 
tstart2=0; 
tend2=data2(2): 
t2=tend2-tstart2; 
("Time window is ’,num2str(t2),’ nanosec’ ] 
("Number of Time Points is ’snum2str(ntp2),”’] 
("Number of subaverages is ’,num2str(nsav2),”] 
back] =zeros(1,ntp2); 
for j=I:ntp2 

back 1(j)=data2(j+4); 
end : 


filename=input(’Enter the ’Background2 filename”’.dat to load :’,’s’); 
eval([’load ’,filename,’.dat’]) 
data3=eval(filename); 
nsav3=data3(1); % Number of Subaverages 
ntp3=data3(4);  % Number of Time Points 
tstart3=0; 
tend3=data3(2): 
t3=tend3-tstart3; 
("Time window is ’,num2str(t3),’ nanosec’] 
["Number of Time Points is ’,num2str(ntp3),”’] 
[Number of subaverages is ’snum2str(nsav3),”’] 
back2=zeros(1,ntp3); 
for j=I:ntp3 

back2(j)=data3(j+4); 
end ; 


% FORMING THE SIGNAL AND NOISE POWER 


ntp=ntp!1; 
s=zeros(1,ntp); 
n=zeros(I,ntp); 
nl=zeros(1,ntp); 


("You will see the signal waveform. Please note the window desired.”’) 
for k=I:ntp 

s(k)=targ(k)-back0(k); 

n(k)=back1(k)-back2(k): 
end; 


% DEFINING TIME AXxIS 
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dt2=10/511; 
t2=0:dt2:10; 


plot(t2,(1000.*(s(1:512)))) 

title (SIGNAL WAVEFORM’); 
xlabel(’Time (nsec)’); 
ylabel(’Measured (mV/m)’); 
grid; 

pause; 

clg; 


x=input(’Enter the begining for time window: °) 
y=input(’Enter the end for time window: °) 


X1=x/0.01953125 ; 

yl=y/0.01953125 ; 

sum 1=0; 

sum2=0; 

for k=xl:y! 
sum l=sum1+((s(k).*2)-(n(k).*2)); 
sum2=sum2+(n(k).*2); 

end; 


% FORMING THE SIGNAL TO NOISE RATIO 


snr=sum 1! /sum2; 
SNR=10.*log10(snr); 


% DEFINING TIME AXIS 


dt=tl/(ntp-1); 
t=0:dt:tend 1; 
dtl=5/255; 
T1=0:dtl:5; 


subplot(211),plot(t2,(1000.*(s(1:512)))) 


ere *,num2str(x),’ to ’,snumd2str(y),’ nsec’)]); 
xlabel(’Time (nsec)’); 

ylabel(*’Measured (mV)’); 

grid; 

pause; 


subplot(212),plot(t2,(1000.*(n(1:512)))) 
title’’/(NOISE WAVEFORM’); 
Xlabel(’Time (nsec)’); 

ylabel(’Measured (mV)’); 

grid 

meta snr; 

pause; 

clg; 
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Of 


» NOISE SPECTRUM 
k=input (‘Enter the smoothing parameter for noise spectrum: ’) 


Meanie, 5: 
(= Gt (Onto -1): 


% UNSMOOTHED NOISE SPECTRUM 


fi —<conn 1): 
Ne=fft(n1,2*ntp); 


plot(f,abs(N(1:ntp))); 

grid; 

xlabel (Frequency (GHz)’); 
ylabel ( Magnitude’); 

title (NOISE SPECTRUM’): 
meta snr; 

pause; 


% TIME SMOOTHED NOISE SPECTRUM 


H=zeros(1,M); 
for 1=1:k 
I(i,:)=n((1-Hi- 1)*M):i*M); 
GG. )=xcorrli ian 
a(iy J=tit(qiti) 2: 
end; 
H=sum(q)/k; 
f1=(12.5/(M/2-1))*(0:(M/2-1)):; 


plot(f1,1000000.*real(H(1:M/2))); 
xlabel (Frequency (GHz)’); 
ylabel (Magnitude * 10*e-6’); 


... ,num2str(k)]); 
grid; 

meta snr; 

pause; 
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APPENDIX D. NOISE PERFORMANCE PROGRAM 
The following Matlab program makes a noise performance 


comparison between the old and new configuration of the TESL. 


% NAME OF THE PROGRAM: NOISE.M 

% WRITTEN BY Lt. ALDO E. BRESANI 

% THIS PROGRAM EVALUATES THE NOISE PERFORMANCE FOR THE 

% TRANSIENT ELECTROMAGNETIC SCATTERING LABORATORY. IT MAKES 
% A COMPARISON BETWEEN OLD AND NEW CONFIGURATION. 

% INPUTS ARE THE CALIBRATION MEASUREMENTS AND THEIR 

% BACKGROUNDS 

clg 

clear 

filename=input(’Enter the ’’New Config.Cal.filename”’.dat to load.... .... aS: ): 


eval([’load ’,filename,’.dat’]) 
datal=eval(filename); 
nsavl=datal(1); % Number of Subaverages 
ntpl =datal(4); % Number of Time Points 
tstart1=0; 
tendl=data!(2); 
tl=tend!-tstartl; 
[Time window is ’,num2str(t1),’ nanosec’] 
[Number of Time Points is ’,snum2str(ntp! ),’’] 
[Number of subaverages is ’,num2str(nsav] ),’’] 
call=zeros(1,ntp!); 
for j=I:ntpl 

call(j)=datal(j+4); 
end ; 
filename=input(’Enter the "Old Config.Cal.filename’’.dat to load.... ....:’,’s’); 
eval([load ’,filename,’.dat’]) 
data2=eval(filename); 
nsav2=data2(1); % Number of Subaverages 
ntp2=data2(4); % Number of Time Points 
tstart2=0; 
tend2=data2(2); 
t2=tend2-tstart2; 
[Time window is ’,num2str(t2),’ nanosec’] 
(’Number of Time Points is ’,snum2str(ntp2),’’] 
("Number of subaverages is ’,num2str(nsav2),’’] 
cal2=zeros(1 ,ntp2); 
for j=I:ntp2 

cal2(;)=data2(j+4); 
end ; 
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filename=input(’Enter the ’’New Config.Back.filename”’.dat to load.... ....:’,’s’); 
eval({’'load ’,filename,’.dat’}) 
data3=eval(filename); 
nsav3=data3(1); % Number of Subaverages 
ntp3=data3(4); % Number of Time Points 
tstart3=0; 
tend3=data3(2); 
t3=tend3-tstart3; 
(Time window ts ’,num2str(t3),’ nanosec’] 
{’Number of Time Points is ’,num2str(ntp3),’’] 
{Number of subaverages is ’,num2str(nsav3),’’] 
back l=zeros(1,ntp3); 
for j=l:ntp3 

back 1(j)=data3(j+4); 
end ; 
filename=input(’Enter the ’*Old Config.Back.filename”.dat to load.... ....:’,’s’); 
eval({’load ’,filename,’.dat’]) 
data4=eval(filename); 
nsav4=data4(1); % Number of Subaverages 
ntp4=data4(4); % Number of Time Points 
tstart4=0; 
tend4=data4(2): 
t4=tend4-tstart4; 
(Time window is ’,num2str(t4),’ nanosec’] 
{’Number of Time Points is ’,num2str(ntp4),”] 
("Number of subaverages is ’,num2str(nsav4),’’] 
back2=zeros(1,ntp4); 
for j=l:ntp4 

back2(j)=data4(j+4); 
end ; 


ntp=notp!; 

sub] =zeros(1,ntp); 

sub2=zeros(1 ,ntp); 

for k=1:1024 
sub1(k)=call(k)-back1(k); 
sub2(k )=cal2(k)-back2(k); 

end; 


sum 1=0; 
for j=513:1024 

sum 1=sum1+{sub1(j))*2; 
end 
x=sum 1/512; 
rmsn=(x)*.5; 
rmsn1=1000*rmsn; 


sum2=0; 


for j=513:1024 
sum2=sum2+sub2(j))*2; 
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end 

y=sum2/512; 
rmso=(y)*.5; 
rmso!1=1000*rmso; 
dt=tl/(ntp-1); 


% DEFINING TIME AXIS 


t=0:dt:tend1; 
dtl=5/255; 
T1=0:dt1:5; 
dt2=10/511; 
t2=10:dt2:20; 


subplot(211),plot(t,(1000.*(sub!))) 

title( NEW CONFIGURATION SUBTRACTED CALIBRATION WAVEFORM’); 
xlabel(’Time (nsec)’); 

ylabel(’7Measured (mV)”); 

grid 

pause; 


subplot(212),plot(t,(1000.*(sub2))) 

title’ OLD CONFIGURATION SUBTRACTED CALIBRATION WAVEFORM’); 
xlabel(’Time (nsec)’); 

ylabel(’ Measured (mV)’); 

grid 

meta noise; 

pause; 

clg 


subplot(211),plot(t2,(1000.*(sub1(513:1024)))) 

title(? NOISE WAVEFORM NEW CONFIGURATION TESL.,...... 
ae RMS VOLTAGE: ’,num2str(rmsn1),’ mV’}); 

xlabel(’Time (nsec)’); 

ylabel(’Measured (mV)”’); 

grid 

pause; 


subplot(212),plot(t2,(1000.*(sub2(513:1024)))) 

title(? NOISE WAVEFORM OLD CONFIGURATION TESL....... 
Se RMS VOLTAGE: ’,num2str(rmso1),’ mV’}); 

xlabel(’Time (nsec)’); 

ylabel(’Measured (mV)’); 

grid 

meta noise; 

pause; 

clg 

end 


ee 
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